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Classical tools of nonlinear dynamics are used to study the highly excited 
vibrations of small molecules. For effective Hamiltonians with one polyad number 
(approximate constant of motion), previously developed methods locate new anhar- 
monic modes using the critical points in the reduced classical phase space. Theoretical 
arguments are given for generalizing the method to more than one polyad number. 
As the simplest classical invariant structure, critical points of the reduced phase space 
are solved analytically without relying on either integrating trajectories or visual in- 
spection. These critical points, especially those that are hnearly stable, are expected 
to indicate regions with the same type of classical dynamics as well as quantum modes 
of vibration. 

The pure bending Hamiltonian of acetylene (C2H2) is analyzed to demonstrate 
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the effectiveness of critical points analysis. Four families of critical points are born in 
distinct bifurcations, each corresponding to a novel anharmonic mode. These modes 
are visualized with custom computer-generated animations. Their origin and nature 
are qualitatively explained through separate consideration of DD-I and I resonance 
alone. Quantum survival probability verifies that the Local and Counter Rotator 
modes are the stable modes of vibration at high excitation. 

The same analysis is extended for the first time to the acetylene stretch-bend 
system, which has never been analyzed classically with all the resonance couplings. 
Prehminary results are obtained for the polyad series containing the C-H stretch 
overtones. The local C-H stretch critical points family, induced by the stretch-stretch 
(Xii/as) resonance, is located and shown to bifurcate into at least 4 new families when 
the stretch-bend resonances are included. The new families indicate that the mixing 
between the stretch and bend may result in novel vibrational modes. 

This dissertation includes my previously co-authored materials. 
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CHAPTER I 

INTRODUCTION 

1.1 Molecular Dynamics Encoded In Spectra 

In its earlier years, spectroscopy was an important tool in determining the 
structure of molecules. Today, the equilibrium structure, as well as the related spec- 
troscopic constants, continues to be actively pursued. However, the accumulated 
theory, techniques and data have also enabled researchers to use spectroscopy as a 
probe of dynamical processes, such as collisions, energy transfer, and chemical reac- 
tions. 

The question is how to decode dynamical information from the experimen- 
tal or calculated spectra. Spectra recorded in either the frequency or time domain 
should reflect the same physical behavior. These domains are formally connected by 
a Fourier transform Electromagnetic radiation affects molecules on at least four 
levels in the order of increasing energy: the nuclear spin, rotational, vibrational and 
electronic Degrees Of Freedom (DOF). The interactions among even these four levels 
within a molecule make it nontrivial to analyze the total dynamics. For chemists, 
the vibrational DOF are of special interest since they are closely related to chemical 
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reaction processes [2]. The focus of this thesis is the vibrational dynamics of small 
polyatomic molecules (3-4 atoms). On the one hand, they are much simpler than 
most organic or biological molecule, where explicit treatment of all DOF is usually 
difficult if not impossible. On the other hand, these small molecules already exhibit 
a rich range of dynamical behavior, especially among the highly excited states, that 
is yet to be fully understood 

1.2 Modes of Vibration 

We start with the frequency-domain spectra. Each resolved level in the spec- 
trum corresponds to one quantum eigenstate of the molecular Hamiltonian. The 
decoding of dynamics involves assignment of these levels with quantum numbers cor- 
responding to their modes of vibration. The word "mode" has been rather liberally 
used to designate patterns in vibrations. The conventional meaning is that there are 
certain coordinates in which the vibration appears particularly simple. Especially, it 
refers to the case when the vibration can be separated into independent periodic os- 
cillations in these coordinates. Such a separation is crucial for comprehending molec- 
ular vibrations that may be too complex for direct visualization. Therefore, mode 
designation is more than mere labeling because it reflects the pattern of underlying 
dynamics. 

A well-known case is the normal modes. When the amplitudes of vibration 
remain small, the force between any two atoms can be approximated as harmonic, i.e. 
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proportional to the displacement from equilibrium configuration. Diagonalization of 
the resulting force constant matrix yields 3A^ — 6 (or 3A^ — 5 for linear molecules) 
normal mode coordinates, N being the number of atoms in the molecule |3]. The 
classical motion of the atoms (determined by Newtonian mechanics) is decomposable 
into oscillations along these normal mode coordinates with characteristic frequencies. 
In the quantum system, the nodes of the eigenstates are distributed along the normal 
coordinates (each intersecting the coordinate perpendicularly). A set of normal mode 
quantum numbers can be assigned from counting the nodes. 

At higher energy, as the vibrational amplitudes increase, inter-atomic forces 
deviate from the harmonic approximation. This requires the inclusion of additional 
terms, such as anharmonic coefficients and resonance couplings, in describing the 
vibrational dynamics. The former distorts the normal modes, while the latter not 
only distort but also mix them. Eventually, the mixing would become so extensive 
that the normal mode picture ceases to be a valid representation of the dynamics. 

It is believed that at sufficiently high energy, the vibrational dynamics enters 
the "bag of atoms" regime which is characterized by a general lack of periodicity. 
In practice, however, this regime may not be easily attained. In at least some and 
perhaps most molecules, even when the normal modes break down dynamics could 
acquire new types of periodicity corresponding to new modes of vibration. For exam- 
ple, in II2O 5] and O3 ^6], the stretching states exhibit a transition from normal to 
local mode behavior as the energy of excitation is increased. In Fig. II. H the normal 
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symmetric and antisymmetric modes have two A-B bonds stretch in phase and out of 
phase by tt, respectively. In the local modes, only one of the A-B bonds is vibrating. 
The more excited states can be assigned in terms of local mode quantum numbers, but 
not the normal mode ones. Here the local mode is not just another way to describe 
the molecule, but a special one that directly reflects the underlying dynamics. More 
exotic modes, such as the "precessional" mode [7j and more complicated collective 
motions |H], are also known to exist. 





FIGURE 1.1 Normal and local stretching modes in an ABA molecule. In the normal 
symmetric (s) and antisymmetric (a) modes, the two A-B bonds vibrate in concert. 
In contrast, the local modes 1 {£i) and 2 (£2) have most vibrational amplitude in only 
one bond. 
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These vibrational modes strongly influence the molecules' chemical behavior. 
When a molecule is between collisions and free of other external interactions, its 
intrinsic reaction rate, if one exists, is determined by the flow of vibrational energy 
within the molecule, a process termed Intramolecular Vibrational Relaxation (IVR) 
|3|9j ■ The complete characterization of IVR requires a "map" of all possible paths 
of energy flow. Yet, in actual modeling of reactions, such a detailed description is 
often simplifled into a statistical treatment. A most popular assumption is the Rice- 
Ramsperger-Kassel-Marcus (RRKM) model, which assumes unrestricted energy flow 
among all vibrational DOF ^Uj- This model is assumed increasingly valid at high 
internal energy and coupling, especially when the corresponding classical dynamics is 
chaotic jTTj. 

In contrast, energy in a stable vibrational mode remains trapped in a small 
fraction of all the possible energy distributions. IVR and reaction probability would 
be very different from the RRKM limit when the molecule is in such a state. Then, 
one has to be cautious about using a statistical model. These stable modes may be of 
great use in coherently controlling chemical reactions [12], because excitation energy 
in these modes may remain long enough for optical manipulation. 

Since the famous Fermi-Ulam-Pasta simulation in 1953 molecular systems 
have served as an application for the mathematical theory of nonlinear dynamics, and 
a motivation for its continued development. The past 50 years have seen an increasing 
interest in applying classical nonlinear mechanics in studying dynamics of microscopic 
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systems. 

1.3 Goal and Structure of Thesis 

The identification and assignment of vibrational modes are crucial for analyz- 
ing the dynamics of highly excited molecules. For over a decade, the Kellman group 
has been using the mathematical tools of bifurcation and catastrophe theory J3] 
to study classical and quantum behavior in the vibration of small molecules. The 
group developed the method of critical points analysis in order to uncover the non- 
linear, anharmonic modes at high excitation via analytical detection (as opposed to 
numerical search). Its earliest formulation focused on the integrable two-oscillator- 
single-resonance systems, according to the critical point (s) in the reduced classical 
phase space [1^]. This was later extended to coupled 3-oscillator systems, which are 
nonintegrable [T7j. These two methods are reviewed in references |18|19j . This thesis 
is a further generalization that takes into account Hamiltonians with multiple polyad 
numbers. The result is a DOF-independent method of generalized critical points anal- 
ysis based on a more rigorous foundation. The vibrational dynamics, especially the 
pure bending subsystem of acetylene (C2H2) is studied using the generalized critical 
points analysis. 

Chapter 2 reviews the relevant background information: the effective quan- 
tum Hamiltonian and polyads, basic tools in classical mechanics, and topics on the 
quantum-classical correspondence. 
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Chapter 3 first describes existing procedures of critical points analysis for 
treating integrable and nonintegrable systems with one polyad number. Then general 
considerations are carried out with regard to the behavior of critical points of an 
arbitrary effective Hamiltonian with polyad number (s). 

Chapter 4 applies this analysis to the C2H2 pure bending system. The re- 
sult reveals the existence of new families of critical points, born in bifurcations at 
increasing energy. The result is compared to those from the separate consideration 
of single resonances, giving a qualitative description of the nature of the new critical 
points. The stable families of critical points correspond to stable quantum modes of 
vibration. 

Chapter 5 extends the analysis to the stretch-bend acetylene system including 
all the resonances. The fate of the C-H normal stretch mode overtones is considered 
under increasing excitation. The preliminary result suggests that a chain of bifurca- 
tions first create a local C-H stretch mode, which then bifurcates into more complex 
stretch-bend modes. 
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Chapter 6 summarizes contributions made in this thesis, and discusses possible 
future directions of research. 



GLOSSARY 



DOF Degrees Of Freedom 

IVR Intramolecular Vibrational Redistribution 

PES Potential Energy Surface 

PO Periodic Orbit (s) 

PPS Polyad Phase Sphere(s) 

RRKM Rice-Ramsperger-Kassel-Marcus model 

SOS Surface(s) Of Section 

ZOS Zero-Order State(s) 



10 



CHAPTER II 

BACKGROUND INFORMATION 

2.1 Effective Quantum Hamiltonian and Polyads 

The molecular vibrational spectrum is often modeled by an effective Hamilto- 
nian obtained from fitting the spectral levels ^. Taking h = 1, the general form of an 
effective Hamiltonian is: 



H.ff = Ho+yv^^ + y:vr+... (2.1) 



with 



H,=J2uJ^U + ^^ + J2 + + (2-2) 

i ^ ^ i,j;i<j ^ ^ ^ 

vt =v^A{dina,nd,r + (a])"(4)"(«on (2.4) 

The Hq term is in the form of a Dunham expansion [21^. fii is a zero-order 
mode (e.g. normal or local modes) number operator, whose eigenvalue is rij. di is the 

^The fit is to either the resolved experimental spectra or the spectra from theo- 
retical calculations l20l. 
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degeneracy of mode i: 1 for non-degenerate modes and 2 for doubly degenerate modes 
(such as the bending of a hnear molecule). Each V2'' term in eqn. (|2.1|) represents 
a resonance that couples the modes i and j. It exchanges m quanta in mode i and 
n quanta in mode j. V^-'^ acts in a similar manner among three modes i, j and k. 
The operators a\, cii have matrix elements identical to those of harmonic raising and 
lowering operators, i.e. 

Vn~Tl\ni + l) (2.5) 
v^l^i-l) (2.6) 
ni\ni) = rii \ni) (2.7) 

|ni,n2,--- ,?T-Ar) comprise a set of eigenstates of -f^o- They are referred to as the 
Zero Order States (ZOS). In the basis spanned by the ZOS, the matrix form of Heff 
is obtained from eqns. ()2.5l2.6p . Diagonalization of this matrix yields quantum 
eigenfunctions in terms of the ZOS. In order to compare these eigenf unctions to the 
molecular coordinate space (such as bond length and angle), one needs to assume for 
each rii an oscillator model, e.g. of harmonic [221 Morse [221 type. 

The resonance coupling terms V cause the ZOS to mix in the eigenfunctions. 
The quantum numbers Ui then are no longer good quantum numbers. However, 
certain linear combinations of them, known as the polyad numbers, may remain con- 
served in the fitting Hamiltonian (and approximately conserved in the exact molecular 
Hamiltonian) . In triatomic molecules like H2O, there is often an approximate 2:1 fre- 



dl\ni) = 
ai\ni) = 
aldilrii) = 
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quency ratio between one normal stretching mode rii and the bending mode n2 
^. This ratio leads to the inclusion of a Fermi resonance term in Heff- The Fermi 
resonance (1) takes one quantum out of rii and adds two quanta to n2, and (2) takes 
two quanta out of 77-2 and adds one quantum to rii. When the equilibrium configura- 
tion of the molecule is non-linear, the bending n2 mode is singly degenerate (^2 = !)• 
The effective two-mode Hamiltonian with Fermi resonance is 

Hpermi =t^l(i^l + ^) + ^2(1^2 + ^) + ^ui^l'l + + Xuini + i)(n2 + ^) 

+ 2:22(^2 + ^? + VFermi[a\ia2y + (0^)^01] (2.8) 

The polyad number 

P = 2ni + n2 (2.9) 
remains conserved since it commutes with the Hamiltonian 

[P, Hpermi] = PHpermi ~ HpgrmiP = (2-10) 

In the quantum Hamiltonian, the presence of P means that the resonance 
coupling only couples ZOS with the same polyad number. For example, there are 
four ZOS \ni,n2) interacting within P = 3: 

|3,0) ^ |2,2) ^ |1,4) ^ |0,6) 



^Meanwhile, the other non-interacting stretching normal mode can be regarded 
as a "spectator". Because the number of quanta in it is constant, it is absorbed into 
the other parameters when we only consider a specific ns manifold. 
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States belong to the same polyad appear in clusters in the spectra, as illustrated in 
Fig. 12.11 Experimentally, it is often the observation of such clustering that leads to 
the adoption of a polyad model 
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|3,0> 
I 2J> 

I M:> 

— I — 1 2.0 > 

— — 

— :: — I o,4> 

— [Q^> 



P=0 P=l P=2 P=3 



FIGURE 2.1 Polyad structures in spectra, a schematic illustration. The states with 
the same polyad number P may strongly interact with each other due to their close- 
ness. Interpolyad couplings, on the other hand, are relatively weak due to the large 
energy spacing between polyads. 
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Under the polyad model there is no couphng between one polyad and another. 
The Hamiltonian matrix therefore can be separated into blocks, each containing the 
intra-polyad couplings. The blocks can be individually diagonalized, which substan- 
tially reduces the amount of computation involved. In the time domain, the existence 
of polyad(s) imposes an approximate restriction on the energy flow: IVR can only 
occur within the same polyad [221 • 

A systematic method to locate polyad numbers was devised by Kellman j27j. 
The method is closely related to an earlier van Vleck perturbation study by Fried 
and Ezra [2H]- Let there be zero-order vibrational modes, excluding the spectator 
ones. Each resonance term V can be represented by a resonance vector in the A^- 
dimensional linear space Vi = {ni,n2, ■ ■ -nN}- AH the resonance vectors Vi in Heff, 
taken together, form a linear subspace with M dimensions (M < A^). Orthogonal to 
this subspace is another {N—M) dimensional subspace from which the {N—M) polyad 
numbers are found. This is graphically illustrated in Fig. 12.21 for A^ = 3, M = 2. The 
coefficients in the polyad number can be taken to be any set of linearly independent 
vectors Pj, which span the (A^ — M) dimensional subspace. In the Fermi resonance 
system, for example, the resonance vector Vpermi = {1, —2} gives the polyad number 
P = 2ni + n2, in accordance with the vector P = {2, 1} orthogonal to Vpermi- 
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FIGURE 2.2 Locating a polyad number from resonance vectors, a schematic illustra- 
tion adapted from Fig. 1 of [TH] . 
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Although the total number of polyad numbers is fixed, each of them is not 
uniquely defined. There is the liberty of multiplying Pj by an arbitrary factor or 
linearly combining any number of Pj. The usual choice is to match their coefficients 
to the approximate integer ratios among the zero-order mode frequencies, as these 
ratios lead to the inclusion of the respective resonance terms in Heff in the first place 

m- 

Even when the more comprehensive Potential Energy Surface (PES) is avail- 
able, effective Hamiltonians are often constructed from the PES with perturbative 
methods 130131 j in order to simplify the subsequent analysis. The effective Hamil- 
tonian is not only a more reliable model (compared to PES) for the highly excited 
states in triatomic or larger molecules, but also easily gives the useful insight of the 
polyad structure. 

The conservation of polyad numbers is never exact ^. The degree of their 
conservation can be estimated by the uncertainty relationship AE ■ At > h. Spectral 
data recorded at low frequency resolution (larger AE) decodes dynamics at shorter 
timescale (smaller At), and vice versa. In molecules, spectral peaks well described 
by a polyad Hamiltonian may break into finer structures when scrutinized at finer 
frequency resolution. This is caused by the small coupling terms not included in the 

^Physical laws of rigorous conservation are based on fundamental symmetry, 
known as Noether's Theorem (see Chapter 12.7 of |n2|)- As an example, the con- 
servation of energy and linear, angular momentum result from the homogeneity of 
time, space and the isotropy of space, respectively. This is not the case with polyad 
numbers. 
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effective Hamiltonian, wliicli exercise tlieir effects (including tfie breaking down of tfie 
polyad number) at longer time scales jHSl- As an example, the acetylene pure bending 
Hamiltonian of Chapter 4 is fitted to spectra recorded with a resolution of 2 cm~^ or 
finer. The corresponding uncertainty At = 2.6 picosecond is much longer than the 
bending vibration period (50 femtoseconds). Since the polyad structure is still present 
at 2.6ps time scale, the polyad numbers predicted from the effective Hamiltonian can 
be assumed valid at the same time scale or longer. 

2.2 Basic Concepts in Classical Mechanics 

Although microscopic systems are governed by quantum mechanics, classical 
mechanics continues to be an important tool in understanding molecular processes 
due to the following fundamental and empirical reasons. (1) Quantum mechanics is 
built upon classical mechanics, as opposed to being a self-consistent theory. Various 
semiclassical methods serve as a bridge between the quantum and classical worlds. 
(2) Even when the classical description is not exact, it provides an intuitive tool for 
the human researcher, whose perceptions are unfortunately macroscopic and there- 
fore classical, to understand the microscopic phenomena. (3) In large and/or highly 
excited systems, treating the whole system quantum mechanically can become chal- 
lenging, making classical and semiclassical methods useful supplements. 

Below we discuss some basic concepts in classical mechanics (of the Hamilto- 
nian formulation) that are pertinent to the topic of this thesis. The reader is referred 
to Tabor j^l] for a general introduction to classical mechanics with emphasis on the 
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nonlinear dynamics. The textbook by Goldstein [32] may serve as a more compre- 
hensive reference. 

2.2.1 Heisenberg's Correspondence Principle 

Heisenberg's Correspondence Principle provides an important connection be- 
tween the quantum and classical worlds . It relates raising and lowering operators 
in quantum mechanics to Fourier components of action-angle variables in classical 
mechanics pUj : 



The A^-mode quantum Hamiltonian Hejf is mapped to an N DOF classical 
Hamiltonian He/f in canonical variables {Ii,(f)i). The hi terms in transform as 



Hpermi = (^ih + ^^2^2 + Xu^l + Xi2hl2 + X22I2 + '^yFermiJ Ifh COs[0i - 202](2.13) 



Classically, the action-angle variables are defined as the conserved action 
and conjugate angle of the zero-order oscillator of mode i. In particular, when 
the oscillator is 1-dimensional and harmonic, they can be related to the Cartesian 
coordinate and momentum by 






(2.12) 



Hence, substitution of eqn. ()2.1H) into ()2.8|) yields 






(2.14) 
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2.2.2 Hamiltonian Classical Dynamics 

The (/j, (pi) variables form a set of canonically conjugate coordinates in Hamil- 
tonian mechanics. Their time evolution (also known as equations of motion) has the 
elegant form 

dlj dH 

0. = ^ = ^ (2.16) 
^ dt dli ^ ' 

Once the initial condition {/^(O), 0j(O)} is given, the subsequent solution {hit), (piit)} 
is determined by integrating eqns. (j2.15l2.16p . {/j(t), 0j(t)} is known as a phase 
space trajectory or trajectory. When the Hamiltonian H is independent of time, a 
trajectory cannot intersect with itself in phase space, although it could retrace the 
same closed orbit repeatedly. 

The equations of motion remain formally invariant. After a transformation 
between two sets of canonical coordinates (called a canonical transformation), e.g. 
(Ji, (pi) —> {Ji, $i), we have: 

2.2.3 Constants of Motion 

In a time-independent Hamiltonian, the energy is conserved, i.e. is a constant 
of motion. Additional constants of motion may be present due to the polyad numbers. 
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In the Fermi resonance system, the polyad number P = 2ni + 712 corresponds to a 
constant of motion through eqn. (j2.12|) : 

I = 2h + h = P + ^ (2.19) 

Like the quantum commutator in eqn. ()2.1Up . the Poisson bracket between I and 
H Fermi also vauishcs. In order for the angle conjugate to / to satisfy 

j = -^El^pm = (2.20) 
0(7 

9 does not appear explicitly in the Hamiltonian. Such a variable is known as a cyclic 
angle [33]. This property leads naturally to a canonical transformation 

(/l,01,/2,02)-(/,^,/.,Vl/) 



with 



^ Jh + h ^ ^=01 + 202 (2.21) 



2 

2/1 - h 



m - 202 (2.22) 



In the new coordinates, the classical Hamiltonian in eqn. ()2.13p is expressed as 

Hperrm =0Ol{I + /.) + CU2(/ - L,) + Xii(/ + hf + Xi2{I^ - ll) + X22(/ - Izf 



+ 2VF,rmi^{I + - h)cOS^ (2.23) 

Since its value does not change with time, / can be regarded as an external parameter. 
9 is absent from the Hamiltonian, and has limited physical significance. The non- 
trivial part of the dynamics is captured in a 2-dimensional phase space (J^, \&), called 
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the reduced phase space. In general, in an DOF Hamiltonian with (A^ — M) 
constants of motion, the phase space can be reduced from 2N to 2M dimensions 
by a similar transformation. The details of such transformations are discussed in 
Appendix A. 

A system is called integrable is there are as many constants of motion as the 
number of DOF. In an integrable system, it is possible to express the Hamiltonian 
in constants of motion only (without their cyclic angles). Then the equations 
of motion can be solved analj^ically without recourse to numerical integration. A 
Hamiltonian with 1 DOF is always integrable when the energy is conserved. 

2.2.4 Invariant Phase Space Structures 

An invariant phase space structure is defined as any lower-dimensional subset 
of the phase space that is mapped onto itself by the equations of motion. These 
structures are the "landmarks" that delineate regions in phase space with different 
kinds of dynamics. The qualitative description of all these regions is called the phase 
portrait [T^ . 

We are primarily interested in non-integrable systems, whose phase spaces are 
at least 4-dimensional and not easy to visualize directly. In order to illustrate the role 
played by these landmarks, the following analogy seems appropriate. Construction of 
the phase portrait of a new dynamical system can be compared to mapping out the 
hydrologic flow on an unknown continent ^. Invariant phase space structures then 

^This comparison certainly should not be taken literally. For example, the au- 
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serve purposes similar to the watersheds, rivers and lakes et cetera in the exploration 
of the continental water system. In the cartoon map of Fig. 12.31 the general flow 
trends on the surface is well characterized by the landmarks in the map, even though 
they do not account for the fate of every single raindrop, and whole regions may be 
left out of the picture (the areas marked with "?"). When there is no prior knowledge 
about the system, finding these landmarks is the first step of a systematic exploration. 



tonomous flow in the phase space has neither "sources" nor "sinks" . 
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FIGURE 2.3 Hydrologic landmarks on an unknown continent, a schematic illustra- 
tion. Regions with different hydrologic dynamics are summarized by a map with the 
watershed (peaks in the center), rivers and lakes. 
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Critical Points The simplest invariant phase space structure is a critical point [37] . 
These are defined as points where the equations of motion (j2.15l I2.16|) vanish: 



Note that the term "critical points" has also been used referring to where the gradient 
of a given function vanishes and the function may not be related to any dynamical 
property. However, in Hamiltonian systems this definition coincides with the one 
defined above in eqn. ()2.24|) . 

Critical points are the simplest invariant structure because they have the lowest 
dimensionality 0, and because they can be exactly solved for as roots of simultaneous 
equations. 

The stability of a critical point intuitively refers to the dynamics of trajecto- 
ries in its neighborhood. Near a stable (unstable) critical point, the dynamics can 
be compared to that near the minimum (maximum) of a classical potential. In the 
illustration of Fig. 12.41 (a), trajectories near a stable critical point are confined to the 
neighborhood, oscillating with small amplitude. In panel (b), an unstable critical 
point behaves locally like a saddle point, with nearby trajectories deviating exponen- 
tially. Mathematically, the stability of a critical point (or other invariant structures) 
has been defined in various ways to suit different purposes. The most important 
ones are Lyapunov, linear and spectral stabilities. Lyapunov stability implies linear 
stability, which in turn implies spectral stability [S^. The linear stability has been 



dH 

d(j)i 



dH 







(2.24) 
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widely used, as it is easy to calculate. It will be the definition used in this thesis to 
characterize critical points; the exact derivation will be discussed in § 3.3.1. 




FIGURE 2.4 Dynamics near linearly stable and unstable critical points in a 1 DOF 
system. Panel (a): trajectories near a linearly stable critical point, (b) trajectories 
near a linearly unstable critical point. 
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Periodic Orbits A Periodic Orbit (PO) is a trajectory that retraces itself with a 
finite period T: 

{Ii{nT),^i{nT)} = {Ii{0),4>i{0)} with n = 1, 2, 3, . . . (2.25) 

A PO is 1- dimensional invariant structure in phase space. Unlike the critical points, 
the only general way to locate a PO is through an iterative numerical search |^ . 
Invariant Tori Another example of invariant phase space structures is the invariant 
torus. In an N DOF integrable system, if Jj are the constants of motion, then their 
conjugate angles 9i evolve at constant frequencies = fj- for any initial condition. 
An example with A^ = 2 is illustrated below in Fig. 12.51 If 6i : 62 happens to be an 
integer ratio p : q (this condition is known as being commensurable), a trajectory will 
close on itself after q periods in direction 9i and p periods in direction 6*2. The surface 
of the torus is therefore covered by a family of PO ^. 

If the frequencies are not commensurable, a trajectory gradually fills the entire 
toroidal surface without closing on itself for any finite time. These quasiperiodic 
trajectories form a set of nesting N-dimensional tori, filling the entire phase space. 



^See the 3D model on the accompanying CD-ROM. 




FIGURE 2.5 Trajectory on an invariant 2-torus. 
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When small perturbations are added to an integrable Hamiltonian, some of 
its invariant tori are destroyed and others become deformed. This is the conclusion 
according to the Kolmogorov-Arnold-Moser theorem . The existence of invariant 
tori in a non-integrable system indicates regions where the local dynamics resembles 
an integrable one. 

In higher dimensions, there are also the normal hyperbolic invariant manifolds 
, which act as impenetrable barriers in the phase space 02] • In an isolated system, 
the {2N — l)-dimensional constant energy shell is also an invariant structure. 

2.2.5 Bifurcations 

A bifurcation generally refers to any qualitative change in the phase portrait, 
as some external control parameters are being varied [14^ . The "qualitative change" is 
typically labeled by the change in the number and/or stability of the invariant phase 
space structures. The "external control parameters" may be either variable physical 
quantities (such as the energy), or the coefficients of the Hamiltonian. 

Fig. l2.6l illustrates the pitchfork bifurcation of critical points in a 1 DOF system 
As the potential V{x) in the Hamiltonian is continuously deformed, suddenly 
the single well (stable critical point) lifts to a barrier (unstable critical point), and 
two additional wells are born. The overall phase portrait changes accordingly, adding 
two zones corresponding to trajectories "trapped" in the two new minima. 




FIGURE 2.6 Pitchfork bifurcation in 1 DOF system and the associated phase space 
change. The right, middle and left panels are before, at and after the bifurcation 
point. 
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In this thesis, quahtative changes in the classical phase space are tracked by 
bifurcations of critical points in the reduced phase space. The parameters in the 
effective Hamiltonian are regarded as given, and the polyad number (s) is the variable 
control parameter. 

2.2.6 Poincare Surface of Section 

The phase space of a 2 DOF Hamiltonian is 4-dimensional and cannot be 
directly graphed like in Fig. 12.61 However, it may be visualized as a series of 2- 
dimensional slices. This technique is called Poincare Surface Of Section (SOS) 

The typical construction of an SOS proceeds as follows. Let the 4 canonical 
coordinates be (/i, /2, 02)- First, an energy value of interest is determined ^, as 
well as a 2-dimensional dividing surface (e.g. by setting 02 at a constant value). An 
ensemble of trajectories at this energy and starting on the dividing surface is then 
calculated. Their intersections with the dividing surface are recorded by two of the 
other independent coordinates (e.g. /i,0i), as illustrated in Fig. 12.71 Due to time- 
reversal symmetry, it is sufficient to record crossings in one direction only, such as by 
letting d(f)2/dt > 0. 



^It can also be some other constant of motion that is held fixed instead of energy- 
see the footnote in 5 3.1.2. 




FIGURE 



2.7 The construction of an SOS. 
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The classical dynamics at this energy is reflected in patterns on the SOS. 
Displayed in Fig. 12.81 are 4 SOS for the Henon-Heiles Hamiltonian i4^, which consists 
of two coupled 1-dimensional oscillators 

Hhh = — ^ + — ^ — + xy - Y (2-26) 

In panel (a) there are two distinct types of trajectories. For each of the red, black, 
blue and green trajectories, the marks remain on a pair of closed curves. It is because 
each of these trajectories lies on an invariant torus, and the two curves on the SOS 
are the result of slicing the torus with a plane (e.g. one that cuts along diameter Ji in 
Fig. 12. 5|) . This kind of quasiperiodic motion is also known as regular. The magenta 
trajectory, on the other hand, randomly fills an area complementary to the regular 
areas, indicating a lack of periodicity. The corresponding random-looking trajectory 
is also known as chaotic. 
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FIGURE 2.8 Regular, mixed and chaotic dynamics from an SOS of the Henon-Heiles 
Hamiltonian. Panel (a) was created in the online program at It displays the 

marks from 5 individual trajectories coded by color. Panels (b)-(d), adapted from 
[34j . are SOS taken at increasing energy. The teardrop-shaped boundary in all panels 
is determined by the conservation of energy. 
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Panels (b)-(d) of Fig. 12.81 are of the same Hamiltonian with increasing energy. 
At the lowest energy (panel b), most of the area on the SOS has a regular pattern. In 
panel (c) there are both regular and chaotic regions. At the highest energy (d), most 
of the SOS is filled with chaotic trajectories. This regular-to-chaotic trend is typical 
for non-integrable systems. In this thesis, we focus on systems with mixed dynamics 
like that depicted in panel (c). The important question here is how to distinguish the 
regular dynamics from of chaos. 

2.3 Quantum- Classical Correspondence 

According to the Bohr correspondence principle, the behavior of a quantum 
system converges to that of the corresponding classical system when Planck's constant 
^ — > 0, or when the quantum number approaches infinity. A more recent theorem 
by Helton and Tabor indicates that in the h —>■ limit, quantum eigenstates must 
localize into phase space regions supporting an "invariant measure", i.e. the invariant 
phase space structures jlH]. 

The correspondence between classical invariant phase space structure and 
quantum wavefunction has also been observed at finite h in a phenomenon called 
localization j46j|. A regular region of classical phase space corresponds to eigenfunc- 
tions with nodes positioned according to the invariant phase space structures. This 
alignment is illustrated in Fig. 12.91 The nodal backbones of wavefunctions closely fol- 
low the PO labeled [r], [B], and [SN]. Wavefunction 4 corresponds to a combination 
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of modes [r] and [B], and the nodes form a rectangular grid distorted along these 
directions. 



38 




FIGURE 2.9 Localization of semiclassical wavef unctions, adapted from [201 ■ The 
curves/lines labeled [B], [SN] and [r] (which coincides with the left edge of upper 
left panel) are the periodic orbits that form the backbones of the wavef unctions. The 
axes are the harmonic normal mode coordinates. 
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In a classically chaotic system, most eigenfunctions lack such regular patterns 
|T7] . However, invariant phase space structures such as PO have been observed to 
have important influence, even when the overall dynamics is chaotic. Examples with 
both stable and unstable PO are known. The latter is called "scarring" by Heller and 
co-workers jlHl- The wavefunctions could be assigned to new anharmonic modes by 
the PO. 

These intriguing observations of quantum-classical correspondence have led to 
renewed interest in using classical mechanics to understand molecular dynamics [49\ . 
The goal is not only finding appropriate phase space structure to explain dynamics in 
an a posteriori manner, but also actively predicting the dynamics from analytically 
detected invariant phase space structures. 
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CHAPTER III 
METHODOLOGY 

In this chapter, first we review two existing schemes of critical points analysis 
for systems with one polyad number in § 3.1, 3.2. Then in § 3.3, additional questions 
concerning arbitrary DOF and multiple polyad numbers are addressed, in order to 
formulate a generalized version of the critical points analysis. The method establishes 
that near a stable critical point in the reduced phase space, classical trajectories are 
quasiperiodic. These critical points therefore indicate the existence of regular modes 
of vibration. 

3.1 Critical Points Analysis of Single m : n Resonance 

This analysis was developed by Kellman et al. for the effective Hamiltonian 
consisting of two zero-order modes coupled by a single resonance. Below, a brief 
overview is given on aspects that will be used in Chapter 4. For a more detailed 
description, the reader is referred to |18p5()j . 

3.1.1 The m : n Resonance Hamiltonian 

In many triatomic molecules, the coupling between two vibrational modes 
1 and 2 (not necessarily normal modes) can be approximated by an m : n type 
resonance. Both m and n are positive integers. Eqn. (j2.8j) is one example with 
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m : n = 1 : 2. For the general m : n resonance the effective Hamiltonian takes the 
following form: 

Hmn = Hoini,n2) + K^n[(a^)™(a2)" + (4)"(«i)'"] (3-1) 

The second term on the right hand side corresponds to a matrix element between ZOS 
\ni,n2) and \ni + m,n2 — n). This coupling destroys both rii, n2 as exact quantum 
numbers, but preserves one polyad number 

rii n2 
m n 

Using eqn. ()2.1ip . a classical Hamiltonian in action-angle variables 
(Ji, 01, J2, ^2) is obtained from Hmn- Then let a be the largest common factor between 
m and n, the following canonical transformation is carried out: 

= 77 --- , ^ = ^^^^ ^ (3.4) 



2 \m 72 ^ 

I is the constant of motion differing from Pmn by a constant, while ^ is a cyclic 
angle. can be regarded as a measure for the extent of mixing between the zero-order 
oscillators 1 and 2, and \& as their relative phase angle. The classical Hamiltonian 
becomes 

Hmn = Ho{I, h) + 2Vmn{I + /.)^(/ - h)^ cos[a^] (3.5) 
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The 2 DOF Hamiltonian of ()3.5p is integrable, since both Hmn and / are 
constants of motion. The classical phase space is reduced to 1 DOF with the equations 
of motion: 



The reduced phase space {h^^) has the same topology as the surface of a 
2-dimensional sphere [Zj. The sphere is called the Polyad Phase Sphere (PPS) ^. As 
shown in Fig. the angle arccos[/z//] is the longitude of the PPS, while \1/ is its 
latitude. The north pole {I^ = 0) and the south pole (Ji = 0) are the mode 1 and 2 
overtones, respectively. At these two points, \& becomes unphysical, since the phase 
angle of an oscillator is ill defined when the action vanishes. 



^Alternatively, (J^, could be regarded as a special SOS in the full phase space 
(/, 9, Iz, Instead of energy, here / is held constant. The dividing surface is defined 
by a constant 6. 



OH 



OH 



(3.6) 



3.1.2 The Polyad Phase Sphere and Critical Points 




FIGURE 3.1 Coordinates on the PPS. The rectangular plane in panel (a) is a Mecartor 
projection of the spherical surface in panel (b). Identical points are labeled by A-D 
in both panels to aid visualization. 
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In an integrable Hamiltonian, each eigenstate of Hmn is associated with an 
invariant torus in (/, 9, 1^, via Einstein-Brilhon-Keller (EBK) quantization [51j. 
On the PPS (which ehminates I and 9), each torus appears as a closed semiclassical 
trajectory. In practice, this trajectory can be well approximated (typically within 1 
cm^^) by simply solving for points on the PPS with the same energy as the quantum 
state 

On the PPS, all semiclassical trajectories are organized by the critical points 
in the reduced phase space: 

dH 



9^ 



(3.7) 



dH , , 

*=gj- = (3.8) 

When there is no resonance (V^n = 0), all trajectories are parallel to the equator 
of the PPS, because Hq has no dependence on \E'. The only critical points are then 
the north and south poles. This corresponds to the trivial case where the eigenstates 
are assigned by the zero-order quantum numbers (ni,ri2)- When the resonance is 
turned on, new critical points may emerge and the old ones may change their sta- 
bilities. The semiclassical trajectories, as well as the quantum states they represent, 
change accordingly. A bifurcation of the critical points signals the birth, death and/or 
transformation of the vibrational modes. 

Fig. 13.21 shows a sample PPS for the HCP molecule (m : n = 1 : 2). Here 
modes 1 and 2 refer to the normal C-P stretch and normal H-C-P bend, respectively. 
In this particular polyad P = rii + n2/2 = 11 there are 12 eigenstates, and their 
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trajectories (labeled 0-11) are evenly spread over the surface of the PPS. The most 
prominent structure on the PPS is a separatrix (dashed line) with the unstable critical 
point [SN] ( "X" ) at the center of its "figure eight" shape. The separatrix is so-named 
because it separates the phase space into three regions: (1) occupied by trajectories 
0-7 which surround the stable critical point [B]; (2) occupied by trajectories 9,11 
which surround the stable critical point [SN]; and (3) occupied by levels 8,10 which 
surround the stable critical point [r] at the north pole. Each level can be assigned 
to two quantum numbers: one is the polyad number P, the other determined by the 
critical point its trajectory surrounds. Since the surface of the PPS is divided, the 
assignment is not uniform for all 12 states. 
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1/2 1 3ro'2 2i 



FIGURE 3.2 PPS and semiclassical trajectories from [201 ■ (a) a view of the PPS. The 
points labeled [r], [5] and [SN] are stable critical points, while [SN] is an unstable 
critical point. The dashed line is the separatrix. Panel (b) presents a cut along the 
great circle defined by \l/ = 0,7r, where 9 G [0, 27r] is a parameter around the great 
circle. Panel (b) also shows the relative energy of eigenstates in this polyad. 
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FIGURE 3.3 Gap in the energy pattern for the eigenstates in Fig 13.21 from j2_QJ . 
Plotted are AE = Ei — E2 vs. Eavg = {Ei + E2)/2 between pairs of levels whose label 
are in the parentheses on the graph. 
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In the full phase space, although /, and \1/ are fixed at the critical points, 
the cychc angle 6 is not. Instead, its value (modulo 2tx) precesses between [0, 27r] at 
a constant frequency. Hence, a critical point in the reduced phase space corresponds 
to a PO in the full phase space. 

3.1.3 Spectral Patterns 

Patterns in the quantum spectra are refiected in the semiclassical trajectories 
on the PPS. First, the ratio between ni and n2 for each state corresponds to the 
time-averaged Iz of its trajectory. In Fig. 13.21 for example, levels 9,11 with Ji -C 
I2 are expected to have a strong bending character. These levels are identified in 
experimental spectra by their large rotational constant [H^ . 

Second, the separatrix on the PPS acts like a barrier in phase space. The 
classical frequency traversing the top of the barrier is expected to drop to zero. The 
quantum equivalence of this frequency is the energy difference Ai? between adjacent 
levels. The pattern of AE therefore exhibits a dip when a separatrix is crossed. 

When there are more than 2 regions on the PPS (e.g. Fig. l3.2j) . Ai? should be 
taken only between eigenstates within the same region on the PPS. If the states are 
sorted by energy alone, levels 9, 11 are intermingled with 8,10. As shown in Fig. 13. 3| 
this choice creates a "zigzag" pattern in panel (a). The smooth dip is recovered in 
panel (b), when the energy differences are taken within the same zone. This resorting 
procedure was first discussed by Svitak et al. in 
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3.1.4 Catastrophe Map 

If Heff includes up to quadratic terms in Hq and Vmn is a constant, all possible 
PPS structures for a given m : n can be further summarized by just two independent 
parameters with the help of catastrophe theory in mathematics ^H] • The PPS up to 
a scaling factor can be reconstructed from these parameters. This 2-parameter space 
(called catastrophe map) is divided into zones for any m : n system [50J, and within 
each zone the PPS have the same qualitative structure. As an example, Fig. 13.41 
displays the catastrophe map and representative PPS for m : n = 1 : 1 ^. In zone 
I, (cases 1, 5, 7, 8 and 9) the spheres share an undivided structure, while in zone II 
(cases 3, 4 and 6) the spheres are each divided by a separatrix. In going from spheres 
1-2-3-4, the bifurcation occurs at sphere 2 where its representative point crosses from 
I (normal mode dynamics) to II (local mode dynamics). 



^3D models of these spheres are also included on the accompanying CD-ROM. 
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m (b) 

FIGURE 3.4 Catastrophe map of 1 : 1 resonance system and associated PPS, adapted 
from ;54; . Panel (a) is the catastrophe map with two independent parameters being 
and A. Panel (b) displays the PPS corresponding to the representative points labeled 
1-9 on panel (a). 
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One limitation of the catastrophe map is that its extension to include high- 
order terms is nontrivial. As shown in § 5.2 of the addition of a single cubic term 
in Hq adds substantial complexity to the catastrophe map. When the high-order terms 
are indeed not ignorable, a simpler alternative using the PPS and spectral patterns 
alone, since they contain the same amount of information. 

3.1.5 Summary 

The steps discussed in § 3.1.1 - 3.1.4 for the single resonance analysis are 
summarized in Fig. 13.51 The 2-dimensional reduced phase space is directly visual- 
ized with the PPS. On the PPS, each quantum state corresponds to a semiclassical 
trajectory. The trajectory can be assigned quantum numbers by the stable critical 
point it surrounds. A separatrix (associated with an unstable critical point) causes a 
"dip" in the energy gap pattern {AE versus Eavg) when the trajectory traverses the 
separatrix. All possible divisions on the PPS with the same m : n resonance can be 
further classified by two parameters on the catastrophe map. 
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FIGURE 3.5 Critical points analysis of the m : n resonance Hamiltonian, adapted 
from Fig. 2 of [HHl with modifications. 
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3.2 Large-Scale Bifurcation Analysis 

In a non-integrable Hamiltonian, the main distinction in the classical phase 
space structure is between the regular and the chaotic regions. Even today, it remains 
a poorly understood field. The most difficult cases have multiple resonances acting 
together, preventing reduction of the dynamics to less than 3 DOF. Lu and Kellman 
proposed the large-scale bifurcation analysis as an extension to these non-integrable 
systems with one polyad number |17|57j . 

The main assumption here is that "The large-scale bifurcation structure is 
defined by the lowest-order periodic orbits and their bifurcations" ^7j. Especially, 
when the reduced phase space has 2 DOF, regions with different types of dynamics can 
be visually recognized on an SOS. Each regular region surrounds a "periodic orbit" 
on the SOS. Here the word "periodic" should not be confused with the continuous 
T in eqn. ()2.25|) for a PO. It refers to the trajectory that appears on the SOS at a 
few discrete points (as opposed to fflling a continuous curve/area). The period is the 
integer number of steps between the returns. Those with period 1 are also known as 
ffxed points on the SOS ^. 

Consider a 3 DOF system with one polyad number, such as the Baggott H2O 
Hamiltonian in |SHI- The polyad number enables one to rewrite the Hamiltonian in 

■^In existing literature, "ffxed points" and "critical points" are often used inter- 
changeably. In this thesis, "ffxed points" is used in the context of a discrete mapping 
(such as an SOS). "Critical points" refer to stationery points in a continuous dynam- 
ical system. 
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a 4- dimensional reduced phase space (/i, V'l, -^2, "^2), plus a conserved action and 
a cyclic angle ips- Dynamics in the reduced phase space can be visualized using a 
series of SOS. Without loss of generality, let the energy and coordinate 'ip2 be held 
constant in the construction of this SOS, and the crossings of trajectories be recorded 
in (Iijipi) space. A fixed point on the resulting SOS has all four action-angle variables 
l2,4'2) constant - therefore it must be a critical point of the reduced phase 

space: 

h = ^i = h = ij2 = (3.9) 

Unless ips has zero frequency, these critical points are closed PO in the full phase 
space. 

These fixed points on the SOS therefore can be found by solving the simul- 
taneous analytic equations ()3.9|) . It avoids numerical integration of many individual 
trajectories, as well as the subsequent problem of classifying their behavior. 

Ref. [T7j solved the bifurcation structure of critical points for several triatomic 
systems. Fig. l3.6l shows the results in H2O. In the limit P — 0, there are 3 branches of 
critical points corresponding to the 3 normal modes in H2O. In the lower right corner 
of Fig. 13.61 the normal bend family is oriented vertically from the origin, while the two 
normal O-H stretch families (on top of each other in this figure) are along the short 
diagonal segment between the origin and point A. As P is increased, resonances cause 
the normal modes to bifurcate (at points A, B, B', etc.) into new families of critical 
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points. These critical points were then used to successfully assign all eigenstates in 
polyad P = 8 to quantum numbers consistent with their vibrational dynamics |IF7] . 




FIGURE 3.6 Large-scale bifurcation structure in H2O, reproduced from Fig. 3 in [17\ . 
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Two subsequent studies |41|59j examined the finer details of tlie pliase space 
structures in tlie same system. Tlieir agreement witli [17j showed the vahdity of using 
critical points to identify the large-scale phase space structures. 

3.3 Generalized Critical Points Analysis 

In § 3.1 and 3.2, the importance of critical points is illustrated for both inte- 
grable and nonintegrable Hamiltonians. Two aspects are worthy of emphasizing: 

(1) The existence of at least one polyad number is crucial for this analysis. In 
reducing the DOF of the classical Hamiltonian, the cyclic angle(s) not explicit in 
the reduced phase space provides time evolution for the critical points in the full 
phase space. In contrast, a critical point defined in the full phase space usually 
conveys little information about the dynamics. For example, although two coupled 
anharmonic oscillators may exhibit a rich range of dynamical behavior, this is not 
apparent from examination of the equilibrium point (usually with no motion in either 
oscillator). With a single polyad number, critical points in the reduced phase space 
are PO in the full phase space. These POs form the "skeletons" of phase space ^j. 
With multiple polyad numbers, the critical points are expected to have the same 
importance, although they now correspond to invariant tori in the full phase space. 

(2) The critical points are found by solving analytically defined equations. Be- 
cause the method does not rely on numerical integration of Hamilton's equations, 
it circumvents the problem induced by unstable/chaotic trajectories. Unlike most 
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existing nonlinear methods, the equations can be extended to arbitrary number of 
DOF without significant change. 

Nevertheless, so far we only considered systems with 1 polyad number and up 
to 3 DOF. The following three points need to be addressed in order to extend the 
analysis to multiple polyad numbers and arbitrary number of DOF. 

1. In a 2 DOF system, the consistency between critical points and large-scale 
phase space structure may be verified by direct inspection, such as through 
SOS. These visual aids are increasingly costly in higher dimensions. Although 
it was suggested that the large-scale bifurcation analysis could be extended to 
> 3 DOF systems with one polyad number ^Hj, it remains unclear how the 
"periodic orbits" defined on a SOS (see § 3.2) can be extended to arbitrary 
DOF. A dimensionality-independent description of the dynamics related to a 
critical point is strongly preferred. 

2. With multiple polyad numbers, the critical points generally have multiple non- 
commensurable frequencies associated with the cyclic angles. Motion at these 
critical points is quasiperiodic in the full phases space, instead of being closed 
PO. To what extent would this difference affect the classical and quantum dy- 
namics of the molecule? 

3. In references |41|57j . the eigenstate assignment was performed through visually 
identifying the localization of the Husimi distribution function of the eigen- 
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states. As both the computation of these semiclassical wavefunctions and the 
visual assignment become impractical in higher dimensions, a more general con- 
sideration of how to assign wavefunction localization is necessary. 

The next three subsections § 3.3.1-3.3.3 discuss these questions in their order. 
The result is a more generalized version of the critical points analysis, which will be 
used in Chapter 4 on the pure bending subsystem of C2H2. 

3.3.1 Reduced Phase Space Trajectory Near Critical Points 

First, we consider an effective Hamiltonian of the most general form. Let the 
Hamiltonian have a total of N modes, M linearly independent resonance vectors, 
and (A^ — M) polyad numbers. The classical Hamiltonian after a suitable canonical 
transformation has 2M action-angle variables spanning the reduced phase space: 

and {N — M) constants of motion and their conjugate cychc angles 

{-Pm+i, ■ ■ ■ , Pn, ^M+i, ■ ■ ■ , ^n} 
Hamilton's equations of motion in the reduced phase space can be written in 



the following matrix form 

/ 



dH 



X 



dXM+l 

dH 

dX2M 

dH 
dxi 

\ / 

\ dXM / 



-Em 
Em 



I dxi * 



dH 

dXM 

dH 

dXM+l 

\ / 

\ dX2M / 



(3.10) 
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with Em being the M x M unit matrix. A critical point Xq in the reduced phase 
space is defined by 2M simultaneous equations: 



dxi) 



] =0 (3.11) 



The linear stability of Xq is defined by the behavior of the linearized equations 
of motion in the neighborhood. Let the point be 

X = Xq + {dxi, . . . , dxi, . . . , dx2M} = Xq + dX (3.12) 

The linearized equations of motion are obtained by expanding dH/dxi on the right 
hand side of eqn. p.lUp into a Taylor series, and keeping only terms linear to the 
displacement 

fdH\ fdH\ V- / d'^H 



- 



Then eqn. flH.lOj) is reduced to the linearized form: 

f -Em \ f d^H 



-X=[ " —— dX = A-dX (3.14) 

dt \ Em J \dxidxjj ^ ^ 

which is a set of homogeneous ordinary differential equations. The standard procedure 

of solving them requires first finding the 2M eigenvalues Aj and eigenvectors Vi of 

matrix A jUU]. The Aj and their respective satisfy 

A-Vi = XiVi (3.15) 

If none of the Ai is zero, the solutions have the following form: 

2M 

X{t) = Xo + ^ e^^* Vi (3.16) 



1=1 
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With Oj being arbitrary complex coefficients. The time evolution of X{t) therefore is 
separable into 2M directions, each indicated by the vector V^. 

The linear stability of Xq is defined in terms of eqn. ()3.14|) . through the 
eigenvalues Aj. In a Hamiltonian system, the conservation of phase space volume 
(Liouville's theorem) leads to the result that Aj always appear in the form of conjugate 
quadruplets (±a ± bi), for which there are four cases described below. 

i. When a pair of A, is purely imaginary (a = 0), all solutions in eqn. (|3.16|) would 
oscillate in the subspace spanned by Vi with a characteristic frequency determined by 
|Aj|. The linear stability in this direction is known as stable, elliptic or (E). 

ii. When a pair of Aj is real (6 = 0), in the Vi subspace all solutions in eqn. ()3.16p 
would be attracted to or repelled from Xq exponentially with time. This direction 
is known as linearly unstable, hyperbolic or (H). The names elliptic and hyperbolic 
originated from the shape of these linearized trajectories (Fig. 12. 4|) . 

iii. When a 7^ 0, 6 7^ 0, the solution contains both oscillating and exponential 
attraction/repulsion components in the subspace spanned by the four Vi correspond- 
ing to the quadruple Aj. In two of the four directions the nearby trajectory "spirals" 
into the critical point, while in the other two directions it "spirals" out of the critical 
point. This stability type is called mixed or (M) [T7] . 

iv. When a pair of Aj = 0, the stability type is degenerate (D). In this case, the 
linearized equations eqn. (j3.14j) become insufficient, and higher-order terms in the 
Taylor expansion are needed to evaluate the stability near a critical point. 
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If all the eigenvalues fall into category i., then the linearized trajectories defined 
by eqn. (|3.16|) oscillate with M distinctive frequencies. Hence, near an all-stable 
critical point, the linearized equations of motion are quasiperiodic. These linearized 
trajectories are expected to resemble the trajectories of the nonlinear Hamiltonian 
Hf,ff for at least a finite time. 

3.3.2 The Presence of Multiple Cyclic Angles 

At a critical point, all the canonical variables are fixed except the (A^ — M) 
cyclic angles 9i. When there is more than one polyad number, the trajectory does 
not close onto itself within a finite time. Otherwise, unless any of their frequen- 
cies becomes zero or commensurable with another, the full phase space trajectory 
is quasiperiodic and restricted to an (A^ — M) dimensional invariant torus. Fig. 13.71 
illustrates the case with A^ = 2, M = 1, which is integrable. The reduced phase 
space (J, is a projection of the full phase space (J, \1', 6). Critical points in (J, \&) 
trace out PO in the full phase space (blue line). A trajectory near the stable critical 
point (green oval at bottom) is a PO in the reduced phase space, and a quasiperiodic 
motion in the full phase space. 
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FIGURE 3.7 Dynamics in the reduced and full phase spaces, a schematic illustration. 
The reduced Hamiltonian is taken to resemble Fig. l2.6l in action-angle variables (J, \E'). 
The angle 9 is the cyclic angle, which evolves (modulo 27r) between [0,27r]. 
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Intuitively, the role these critical points play in the phase space should not 
change whether there are one or more cyclic angles. As an example, consider the case 
of HCP where only 2 of the 3 normal modes are coupled by a Fermi resonance pUj . 
Excitation in the spectator mode 3 (C-H stretching) can be treated as a parameter 
in the effective Hamiltonian. So strictly speaking, there are two polyad numbers and 
cyclic angles: 

P,=m + ^, ^, = 201 + 02 (3.17) 

P2 = ns, 62 = h (3.18) 

Yet, in assigning e.g. the na = 1 states, one could simply use the critical points found 
in this manifold, in spite of the fact that the frequency of 62 is not zero in these 
polyads. 

We argue that the cyclic angles 6i in general represent a trivial aspect of the 
dynamics. As far as quantum assignment is concerned, this is evident if one considers 
the semiclassical quantization procedure. There 6i appear only in a pre-factor with 
the form Ile*^'^' in the resulting wavefunctions j8j. In classical mechanics, also note 
that 6i do not have physical meaning on their own, since the polyad number Pj are 
not uniquely defined (Appendix A). 

Critical points in a reduced phase space, especially those with non-zero fre- 
quencies in the cyclic coordinate, have been known as relative equilibria in mathe- 
matical literature |6H62j . Near a relative equilibrium, classical dynamics in the full 
phase space can be separated into two parts: the group orbit, which is motion along 
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the cyclic angles, and motion in the reduced phase space The latter is a mul- 
tidimensional "slice" transverse to the group orbit [U^. The slice contains all the 
"essential dynamics" in the sense that the full dynamics can be reconstructed 
from a point on the slice and appropriate initial conditions. This provides a further 
argument against making a distinction between systems with one and multiple polyad 
numbers. 

In the field of chemistry, relative equilibria theory has been used to classify 
rotationally excited molecular spectra |65l66j . The total angular momentum J plays 
the same role as the polyad numbers in this thesis. A stable relative equilibrium cor- 
responds to the molecule rotating with a fixed shape. Vibrational modes are defined 
by the normal form of the Hamiltonian in the neighborhood. As J is increased, the 
bifurcations of relative equilibria correspond to predictions of the (as yet unobserved) 
rovibrational spectral patterns. 

3.3.3 Semiclassical Localization Near Critical Points 

An eigenstate may be assigned meaningful quantum numbers based on the 
critical point, if its representation in the same Jj) space is localized near the 
critical point with a well-ordered nodal pattern. The semiclassical wavefunctions 
can be obtained through either phase space representations (e.g. the Wigner or 
Husimi function ^67j), or the space semiclassical quantization proposed by Voth 
and Marcus in [HH]. Examples from both methods are illustrated in Fig. 13.81 In 
panels (a) and (b), the localization occurs around the normal mode critical points at 
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(rii —n2 — 3.5, ■01 = f and y^). In panels (c) and (d), the localization is around the 
critical points at (V'a = ±7r, = ±7r). 
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FIGURE 3.8 Semiclassical localization near critical points. Panels (a), (b) are repro- 
duced from Fig. 3 of jHZI, displaying the two projections of the Husimi function of 
the same eigenstate in action {ni,n2) and angle '?/'2) space, respectively. Panels 
(c) and (d) are reproduced from Fig. 5 of [69j, which display the angle-space repre- 
sentation of two different wavefunctions both localized around {%lJa = ivr,?/';, = ±7r). 
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Consider a local minimum or maximum (together referred to as extremum) in 
the reduced phase space Jj). This extremum point is necessarily a critical point. 
Then if there is a quantum eigenstate whose energy is nearby, intuitively one expects 
the semiclassical representation of the eigenstate in either (^'j, Jj) or (\E'j) space to 
localize near the critical point, simply because of the limited volume of accessible 
phase space into which it can expand. This is illustrated schematically in Fig. 13.91 in 
the case of a minimum. 




FIGURE 3.9 Localization near a minimum in the reduced phase space, which is also 
an elliptic critical point of the Hamiltonian. An eigenstate whose energy is close to 
the minimum must be localized in the nearby phase space. 
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There is no apparent reason why the same argument should not be vahd for 
all choices of semiclassical representation as well as for arbitrary DOF, except for the 
following two scenarios. The localization may be disrupted by quantum tunneling 
when the local extremum is not sufficiently prominent, or there are other local extrema 
nearby with similar energy. 

In § 3.3.2, it was shown that the linearized motion near an all-stable criti- 
cal point is quasiperiodic. If (1) it is a good approximation for the real classical 
trajectories in this region and (2) the region is large enough to support one quantiz- 
ing invariant torus, then the semiclassical wavefunctions may be localized around the 
torus, with M quantum numbers assigned by EBK quantization. Therefore, one could 
expect the all-stable critical points to correspond to quantum modes around which 
semiclassical wavefunctions localize. Other critical points with partial linear stability 
(while unstable in some directions) may also become the center of localization under 
favorable circumstances. 

3.3.4 Summary 

Prom § 3.3.1 - 3.3.3, we can draw the following conclusions about critical points 
in the reduced phase space: 

1. Along the stable directions of a critical point, the linearized classical trajectories 
nearby are quasiperiodic. 
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2. The presence of multiple cyclic angles is not expected to affect the essential part 
of the classical dynamics or semiclassical localization. 

3. If a critical point is also a local extremum in the reduced phase space of i^e//, 
then it is expected to be a center of localization for semiclassical eigenfunctions. 

Therefore, the critical points can be used to assign vibrational modes to the 
quantum spectra. The change in their number and/or stability should correspond to 
the change in birth, death and transformations of the vibrational modes. 
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CHAPTER IV 

BIFURCATION ANALYSIS OF C2H2 BENDS 
4.1 Introduction 

Acetylene (C2H2) is among the most-studied polyatomic molecules in spec- 
troscopy. Its normal mode constants on the ground electronic state (So) have been 
refined over a long period of time jTUI. At increased energy, however, its vibrational 
dynamics is rather complex. Recently the highly excited vibrational states acety- 
lene became more accessible due to new techniques such as Stimulated Emission 
Pumping. The observations include additional spectral features under enhanced res- 
olution [32], signatures of both regularity and chaos in level statistics ^Tlj, and the 
emergence of local modes in H-C-C bending |3U|ti9l72l73l74|75] and C-H stretching 
dynamics |26|76|77] . The stretch-bend system has also been investigated in a few 
studies |78|79|8(J|8H82|83|84j . Theoretical analysis of these results, however, remains 
challenging: while the coupled vibrational DOF are too numerous for many analytical 
methods, they are not enough to warrant a statistical approach. 

The dynamics of acetylene vibration has an important role in the combustion 
processes. Its interconversion with vinylidene, a marginally stable isomer, affects 
the outcome of reaction rate modeling jHH] . Fig. 14.11 illustrates the geometry and 
energy changes during the isomerization process. The transition state is planar and 
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involves mostly one C-H bond in acetylene bending over. The highest bending levels 
recorded by Field et al. are believed to be within 5,000 cm~^ of the transition state 
[Sni • Decoding the dynamics hidden in these spectra is expected to shed light on this 
isomerization process. 
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FIGURE 4.1 Acetylene- vinylidene isomerization on the Sq electronic surface. The 
energy and configurations are averaged over compiled ab initio results in Table 1 of 
ISIj. Note that 1 eV = 8080 cm-\ 
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For modeling highly excited states (~ 10,000 cm~^), the effective Hamiltonians 
from fitting spectra are more reliable than any current ab initio PES. The best-known 
PES ^88j is only qualitative in reproducing the experimental spectra in this energy 
range [HH] . A refined PES recently published by Bowman et al. is still not as accurate 
direct fit to the spectra |90|91j . 

Near its linear equilibrium configuration, acetylene has 3A^ — 5 = 7 normal 
modes as shown in Fig. 14.21 Both 1/4, 1/5 bends are doubly degenerate. Two additional 
quantum numbers £4, £5 are used to label the respective vibrational angular momenta. 



fi Sjirmietric 
C-H Stretoli 



Vi) Trans-Bend 



'Vs Antisymmetric 
C-H Stretch 



o — €^0- 



ii, k Vibrational 
Aiigular Kfomenta 



FIGURE 4.2 Normal vibrational modes of C2H2. 
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All existing stretch-bend effective Hamiltonians conserve three polyad num- 
bers: 

Nt = 5ni + 3^2 + 5ra3 -f 714 -f- 715 (4.1a) 
Ns = ni+n2 + ns (4.1b) 
£ = 4 + 4 (4.1c) 

Ni is the "total" quantum number, representing the approximate integer ratio among 
the normal mode frequencies. Using the frequencies of [^2] in units of cm~^, the ratio 
is 

cji : W2 : CJ3 : CJ4 : = 3371.66 : 1974.76 : 3288.75 : 608.50 : 729.18 ^5:3:5:1:1 

Ns is the total number of stretching quanta. In the absence of rotational excitation 

(J = 0), £ denotes the total angular momentum of the molecule, which is always 

conserved. The inversion symmetry (gerade/ungerade) and parity (+/-) are also 

conserved jHSl- Each polyad can therefore be labeled as [Nt, Ng, i]^/^. 

Incidentally, all three polyad numbers remain good for the isotopomer ^^C2H2 

|M] . The doubly deuterated C2D2 conserves A^^ and i [95] . Experimental spectra also 

suggest the existence of polyad structure in the monodeuterated C2HD pHj . 

4.2 C2H2 Pure Bending System 
4.2.1 Quantum Effective Hamiltonian 

Among the [Nt, Ng, i]^^"^ polyads, the pure bending states with {Ng = 0) form 
a separate subsystem. In the remainder of this chapter, these polyads are labeled by 
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[A^6,£]s/" with 

Nk = Nt = n^ + n, (4.2) 

The latest pure bending effective Hamiltonian was produced by Field et al. 
The data used in their fit include both energy and intensity information from FTIR 
and Dispersed Fluorescence spectra [22] ■ The high- lying levels were obtained from 
DF spectra with a frequency resolution of 2 cm~^, and unresolved i = 0/2 states. 
These levels are as much as 15,000 cm~^ above the ground vibrational state. The 
parameters in this Hamiltonian are listed in Table 14.11 

TABLE 4.1 Parameters in pure bending effective Hamiltonian, from [22] ■ The pa- 
rameters are in units of cm~^. 



004 


608.657 


2/555 


0.00955 


UJ5 


729.137 


fl'44 


0.677 


X44 


3.483 


5'45 


6.670 


3^45 


-2.256 


fi'55 


3.535 


X55 


-2.389 


5'45 


-8.574 


1/444 


-0.03060 


^0 
'45 


-6.193 


2/445 


0.0242 


^"445 


0.0304 


2/455 


0.0072 


^545 


0.0110 



Using the normal mode basis \n^,n^^) {rii > as the ZOS, the fitting 
Hamiltonian has a diagonal part Hq and three resonance couplings. 

Hbend = Hq + Vddi + Vi + Vddii (4.3) 
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Ho = uj^nA + ws^-s + Xi^nl + x^^nin^ + x^^nl + Vi^nl + 2/445^4^^5 + Vi^^n^nl 

+ 2/555^5 + ^44^4 + 9idik + g^dl (4.4) 

1. A Darling-Dennison |23 resonance (DD-I): 

{n[\n'i\VDDl\{ni - 2f\{n^ + 2)^^) = ^ [K - ll){n^ + 4 + 2){n^ - 4 + 2)] (4.5) 

2. An ^-resonance: 

(n^/, n^55|F,|ni*^2^ n'^^^) = ^ [{n, ^ h) {n^ ± ^4 + 2) (ng ± 4) (^5 T 4 + 2)]^/' (4.6) 
with i?45 = r^s + r445(n4 - 1) + r545(n5 - 1). 

3. Another Darhng-Dennison resonance (DD-II), with matrix elements smaller than 
those of the previous two: 

lo 

[{n, ± Q{n, ± £4 - 2)(n5 ± 4 + 2)(n5 ± i, + (4.7) 

Fig. 14.31 illustrates the manner these resonances act within a polyad. The 
resonances are shown as lines connecting pairs of ZOS. Vddi couples within each 
column (same ^4,^5), while Vg couples within each row (same 71,4, ns). Vddii in eqn. 
()4.7|) contains couplings along the diagonal directions. With both Vddi and Ve, or 
Vddii alone, all ZOS in a polyad [Nf,, Pfl"^ are connected into an inseparable network. 
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FIGURE 4.3 Resonance couplings within a pure bending polyad [8,0]. Panel (a) 
displays the g states and panel (b) the u states. 
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Alternatively, Hbend can be written with raising/lowering operators, which act 
on the normal mode trans- and cis- 2-dimensional oscillators. The symmetry-adapted 
operators a\^, a\g, a^d, o^ig and d^^, a\g, d^d, d^g ^ are defined according to ^98j as 

aLl^f ) = y^^^f^K-. + 1)^^^^), al>f ) = y^^^f^K-. + It-') (4.8) 

) = ^|^^\{n, - lY^-'). a,,|nf ) = ^J^^\{n, - (4.9) 
The zero-order part of i/;,end can be expressed with 

Ui = dl^did + dlgdig, ii = d\^did - dJ^Sig (4-10) 
and the resonance terms as 

yODI = 'S'45 [0,4^04^05^053 + d4^dO'ig<^\d'Ag] (4-11) 

Vi = R45[(^Md.\gdl^a5g + ai^ciigdl^d^g] (4-12) 

VdDII — ^ [a4rf04d^5dCt5(i + O'AgO'ig'^^g'^^g 

+ duaAgdl^dl^ + d4<;a4gaJ3ay (4.13) 
4.2.2 Classical Hamiltonian 

Next, a classical Hamiltonian is obtained from Hbend using Heisenberg's Corre- 
spondence Principle of eqn. (12.111) . The result is expressed in canonical action-angle 

^The g (gauche/left) and d (droit/right) subscripts are named according to the 
manner they affect the vibrational angular momentum ^i. 
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type variables {i = 4:d,4:g,5d,5g). The actions are related to the zero-order 

quantum numbers by 

hd ={714 + 1 + Q/2, hg = (n4 + 1 - £4)72 (4.14) 

ha = K + 1 + 4)/2, hg = (n5 + 1 - 4)/2 (4.15) 

Corresponding to A*";, and there are two conserved classical actions: 

hd + hg + hd + hg = ni + n^ + '2 = Nb + 2 (4.16) 

hd-hg + hd-hg = ^i + ^^ = ^ (4.17) 

In order to reduce the dimensionality of the phase space, the following canon- 
ical transformation is useful ^: 



hd 


+ hg 


+ hd 


+ hg 


Nk + 2 






2 




2 ' 


hd 


- hg 


+ hd 


- hg 
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~ 2' 


hd 


+ hg 


— hd 


- hg 


n4 — ns 






2 
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hd 


- hg 


— hd 


+ hg 





j.4:d ^Ag I ^bd ^bg ^ n 

K, = = -, 9, 

_ hd + hg — hd — hg _ — Ub 

~ 2 ~ 2 ' 

^ hd - hg - hd + hg ^ ^ ' (4.18 



04d + 0. 


4g + ( 


hd + (j)5g 




2 






45 + < 


hd - 4>5g 




2 




4>4d + 4>. 


43 - ' 


hd — 05g 




2 




04d - 0. 


4g - < 


hd + 4>5g 



^ These new coordinates are defined in the same way as Jacobson et al. in 
The only difference is that our actions are expanded by a factor of 2, while the angles 
are reduced by a factor of 1/2. Such a difference is trivial, except ours are more 
compatible with the single-resonance analysis of § 4.4.1. 
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The resulting classical Hamiltonian is 



H^^^^iKa, Kb, Ja, Jb, 1pa,A) = ^i{Ka + Ja " 1) + ^b{Ka - Ja " 1) 

+ Xii{Ka + Ja- if + X45iKa + Ja - l)(i^a " Ja " 1) + 3^55 (i^a " " 1)^ 
+ vMKa + Ja- if + yU5iKa + Ja - lf{Ka - Ja - 1) 
+ VAb^iKa + Ja- l){Ka - Ja - if + y^5b{Ka - Ja - if 
+ gu{Kb + Jbf + 545(i^fe^ - Jfe^) + g55{Kb - Jbf 



+ 
+ 



2 
2 

1 



{K', - Kt) + {J'a ~ Jtf - mt + KtfJ'a + Jt) - SKaKbJaJb COs[2V'a 



_ Kif + {Jl - Jif - 2{Kl + Kl)(Jl + J2) _ SK^KbJaJb ' cos[2^,] 



+ ^[Rl^ + 2545]{[(Aa + Kbf - [Ja " Jfe)'] COs[2(Va " ^b)] + 
[{Ka - Kbf - {Ja + Jbf] COs[2(V'a + V'fe)]} 



(4.19) 



with 



^Is = ^^45 + '^445(^a + Ja " 2) + r^^^^Ka - Ja ' 2) 



(4.20) 



In order for all the actions l4d, Jaq, Ibd and I^g in eqn. ()4.18|1 to be non-negative, 
additional restriction are applied to the values of Ja, Jb'- 



Ka- Kb>\Ja + Jbl Ka + Kb>\Ja- Jb\ (4.21) 

Since Ka and Kb are constants of motion, their conjugate angles da-, Ob are the 
cychc variables and therefore absent from the Hamiltonian. However, they do evolve 
in time with 



(4.22) 
(4.23) 
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The reduced phase space {Jayipa, Jhii^b) is four- dimensional. The equations of 
motion in it are: 

) dH -OH ,^ 

= ^- ^' = - 1? (^-25) 

The parameters in Table apply to £ = 0,2, which from eqns. ()4.18|) cor- 
respond to f^;, = 0, |, respectively. This chapter will focus primarily on these cases. 
For £ = = 0, eqn. ()4.19|) becomes 

Hbend{Ka, Ja, Jb, ^a, ^h) = Hq + VdDI + Vi + VdDII + V'oDII (4-26) 

with 

= Ui{Ka + J„ - 1) + U^iKa -Ja-l)+ Xu{Ka + Ja ' 1)' 

+ X45{Ka + Ja- l){Ka - Ja - 1) + Xgsli^a ' Ja - if + 2/444(i^a + Ja - if 

+ 2/445(^a + Ja- l)'(i^a " " 1) + l/455(^a + Ja ' l)(i^a - Ja - 1)' 

+ l/555(^a -Ja- if + {QaA " 9^5 + ^755)^' (4-27) 

and 



Vddi — 


^[K + {Jl-Jl) 


''-2Kl{Jl + Jl)]"' 


COs[2lpa\ 


(4.28) 


Vi = 


^[K + {Jl-Jl 


?-2Kl{Jl + Jl)]'/' 


cos[2'?/^fc] 


(4.29) 


Vddii = 




{Ja - Jbf] C0S[2(^, - 




(4.30) 


^DDII = 


hRl^ + 2g,,)[Kl- 


{Ja + Jbf] C0S[2(V^, + 


-A)] 


(4.31) 
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4.3 Critical Points Analysis 



Following the method outlined in Chapter 3, we now explicitly solve for the 
critical points for ^ = in eqn. (|4.26p and their bifurcations with variation of A*";,. 
Four new families of critical points are found at increasing Nb when the initially stable 
normal modes become unstable in distinct bifurcations. Two of the new families, 
namely the Local and Counter Rotator critical points, are linearly bi-stable (EE). 
They correspond to new stable modes of bending vibration. 

Readers not interested in details of the calculation may skip the next subsection 
and go directly to § 4.3.2 for the results. 



The critical points of eqn. (|4.26|) are defined by four simultaneous equations: 



4.3.1 Computational Details 
















(4.32) 



(4.33) 



(4.34) 



a 








(4.35) 



To simplify the notation we let 



A = IK 



'^Ki{Ji+Ji) + {Ji-Jini' 



(4.36) 
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The partial derivatives in eqns. 



become 



= - A545 sinpV',] - -^{Rl^ + 2g^^){[Kl - {Ja - J^f] sm[2(^, - ^6)] + 

[Kl - {Ja + Jb?] sm[2(Va + m} (4.37) 
BH 1 

— ^ = - KRl^ sin[2V'6] + -(iJis + 2545){ [i^a " {Ja - Jb?] sinpCV-a - ^b)\- 

OiPb 4 

_ (j^ + j,)2] sin[2(V., + ^,)]} (4.38) 

+ 2^45 



{{Ja - Jb) COs[2(^a - Vb)] + (Ja + Jb) COs[2(V'a + ^fe)]} 



+ 



r445 - r545 



{2J« Jf, sin[2V^„] sinpV^b] + (K^ _ j2 _ j2) cos[2^„] cos[2V^b]} 



(4.39) 



dH, 



bend 



Jb 



dJb 



2(544 - 545 + §55) Jb - ^{Ka + Ja - Jb){SAb COs[2'0a] + ^45 COS^V-b]) 



A 



-^45 + 2^45 



{{Ja - Jb) COs[2{lPa - V'b)] - {Ja + Jb) COs[2(Va + i^b)]] (4.40) 



A visual inspection of the first two equations ()4.37l 14. 38^ reveals a sufficient 
condition ^ for both of them to vanish: 

sm[2^a] = sinpV^b] = (4.41) 

or 

TTITT TITT 

(V'a,^6) = ( — ,^) withm,n = 0,l,2,3 (4.42) 

Eqn. ()4.42p leads to the 16 combinations of {tpa, '4'b) in Table 14.21 The re- 
maining two equations ()4.39| I4.40|) are then solved for Ja and Jb, with {ifja, ^Jb) held 

^An additional root search was carried out, and no extra critical points were found 
beyond the ones discussed here. 
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at these discrete values. Also, note that eqns. ()4.39t I4.4U|) are invariant with regard 
to one or both the following two transformations 

V'a ^ V'a + TT Ipb ^ i^b + (4.43) 

It is sufficient to consider only one out of the four sets of values on each row in 
Table Here we use: 

(^a,^5) = (O,O),(o,^),g,o),(^,0 (4.44) 
TABLE 4.2 {ipa, ipb) values of bending critical points. 



(^a, 'ipb) 








(0, 0) 


(0, vr) 


(vr, 0) 


(tt, vr) 


(0, f) 


(0, f) 






(f,o) 




(1,0) 




/ TT ■K\ 


/ n 3tt\ 




/ 37r 37r \ 
y 9 ^ 9 1 



The simultaneous eqns. ()4.391 14.40|) may be visually solved using the following 
"pseudopotential" approach. Here it is generalized from the 1-dimensional case in 
[20] to 2-dimensional. When i^bend is viewed as a function of Ja, Jh with the discrete 
i'^ay'^b) values in eqn. ()4.44j) . the solutions satisfying both ()4.39l W^^^ are the "flat 
spots" (e.g. minima, maxima and saddle points), and can be identified visually on 
plots of Hbend (the pseudopotcutial) . Fig. l4.4l depicts the pseudopotentials for polyads 
[6, 0] and [22, 0]. The Ja and Jb values (already constrained by eqn. I4.21|) are divided 
over Ka to scale their range to [—1,1]. In panels (a)-(d), there are no fiat spots in the 
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interior of the {Ja/Ka, Jb/Ki,) space when Ni, — 6. In panels (e)-(h), under each of the 
four (-0(1, ■^b) conditions there is a flat spot (white dot with black rim) corresponding 
to a critical point. This indicates the existence of at least four critical points when 
TV, = 22. 
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However, when Ka is treated as a continuously varied parameter, it is more 
straightforward to solve (|4.391 l4.4Up analytically, instead of locating them graphically 
through the pseudopotentials. Given the (■i/'a, ilJb) values in eqn. ()4.44|) . eqn. ()4.4Up is 
further simplified to: 

-jj^ = M2i9u - 94, + 9,,) - " ^° ' (^45 cos[2^J + R*,, cos[2^,]) 

RI5 + 2(745 cos[2z^,] cos[2z^b] } = (4.45) 



All terms on the right-hand side are proportional to Jb. The other multiplier (the sum 
of terms between the curly brackets) can be numerically shown to be always positive. 
Hence eqn. ()4.35|) has the trivial solution: 

Jt = (4.46) 

The last equation ()4.34j) is then solved analytically on a computer using the software 
Mathematica |99j . 

At each critical point, the linear stability is determined by the eigenvalues of 
the stability matrix described in § 3.3.1. Here the matrix is 4 x 4 in size, and the 
possible stability types are bi-stable (EE), bi- unstable (HH), stable-unstable (EH), 
mixed (MM) and degenerate (D). 

Special consideration of the critical points is required where \Ja\ + \Jb\ = Ka- 
At these points, the values of ipa and ipi, are indeterminate as the denominator A in 
eqns. ()4. 39(100)) vanishes. The ( Ja, ipa, Jb, i^b) coordinate system becomes singular at 
these locations. An alternative coordinate system is required to evaluate (1) whether 




(e) (f) (g) (h) 




FIGURE 4.4 Pseudopotentials of [6,0] and [22,0], as contour plots of -ffbend- The 
gray scale shading is the lightest where Hbend is at the largest. Panels (a-d) depict the 
pseudopotentials 'm.Ni, = Q with the four values of (^a, ip},) in eqn. ()4.44|) . respectively. 
Panels (e), (f), (g) and (h) depict the pseudopotentials in = 22 with the same 
order in (^a, V^fe)- 
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a point is critical point or not, and if yes, (2) its linear stability. The technical 
details are discussed in Appendix B. Only the following four are critical points with 

\Ja\ + \ Jb\ = Ka- 

Ja = ±Ka, Jb = with (EE) stability (4.47) 
J„ = 0, Jb = ±Ka with (MM) stability (4.48) 

4.3.2 Results of the [Nb, 0] Polyads 

Normal Modes at Low Up to Nb = 6, there are two families of critical points 
at Ja = Ka and Ja = —Ka with (EE) stability. They correspond to the normal 
Trans and Cis modes, respectively. This is consistent with the usual assumption that 
vibrations near equilibrium are dominated by the normal modes. 
Bifurcation at Higher Nb With increasing polyad number Nb (as well as increasing 
energy), four bifurcations occur to the Trans and Cis critical points. Four new families 
of critical points called the Local (L), Precessional (Pre), Orthogonal (Orth) and 
Counter Rotator (CR) emerge out of the normal mode families at the points of the 
bifurcations. The analytical solutions of these families are listed in Table 14.31 and 
eqns. fj4.49ll4.to|l . 
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TABLE 4.3 New critical points in [Ni,, 0] polyads. 





Name 




Stability 


L 


Local 


(0,0) 


(EE) 


Orth 


Orthogonal 




(EH) 


Pre 


Precessional 




(HH) 


CR 


Counter Rotator 


/ ir 7r \ 
V2' 2) 


(EE) 



L: Ja = 21.166 - 0.584Ka - 0.2681^4091.902 + 42.376/s:„ - (4.49) 

Orth: Ja = 24.876 - 0.496ii'„ - 0.8227^531.865 + 32.598Ka - (4.50) 

Pre: J„ = 10.920 - 0.564X„ - 0.4603^-276.741 + 93.490^^ - (4.51) 

CR: Ja = -18.312 - 0.536Ka + 0.6344^309.356 + 143.457X„ - (4.52) 

At each of the four bifurcation points, either Trans or Cis changes its stabihty, and 
gives birth to one new family of critical points in the following manner: 

Trans (EE) ^^^zL!^ Trans (EH) + L (EE) (4.53a) 

Cis (EE) Cis (EH) + CR (EE) (4.53b) 

Trans (EH) ^"^^'^^ Trans (HH) + Orth (EH) (4.53c) 

Trans (HH) 3^}^ Trans (EH) + Pre (HH) (4.53d) 

The Trans family undergoes three consecutive bifurcations while the Cis family 
undergoes one. All six families involved in these bifurcations have Jf, — 0, and thus 
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^4 = ^5 = 0. The ratio (Ja/Ka) is presented in Fig. 14.51 as a function of Nj, ^. At each 
bifurcation point, the new family appears at the respective Ja = ^Ka, then migrates 
towards Ja/Ka = with increasing Njy. This trend reflects increasingly equal mixing 
between the Trans and Cis in the new critical points, which will be explained in detail 
below. 



"^In Fig. 14. 5t the bifurcation is calculated up to Nb = 30, although the effective 
Hamiltonian is based on spectra up to only 22. 
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FIGURE 4.5 Critical points and their bifurcations in [Nf,, 0] polyads. The four new 
famihes of critical points are shown on the figure, while the Trans and Cis families 
coincide with the upper and lower edges of the figure, respectively. The stability of 
critical points is indicated on the figure. 
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Trajectories In the Full Phase Space The critical points in the reduced phase 
space have two cyclic angles {9a, 9b) that are absent from Hbend- corresponds 
to an overall translation in time, while 9b corresponds to a rotation of the system 
around the C-C axis jHS]- According to § 3.3.2, in the full phase space such a critical 
point is a quasiperiodic trajectory on a 2-dimensional invariant torus. 

However, the [A*";,, 0] polyads have Kb = and thus the frequency of 9b is 
no longer physical. A proper analogy is that when one "shrinks" to zero the tube 
thickness of the invariant torus, quasiperiodic trajectories on its surface are reduced 
to periodic orbits with a single frequency 9a- At the critical points the frequency 9b 
is zero when we substitute in Kb = Jb = 0: 

■ [dHbend\ ^j, , , X Si^JbCOs[2i)a]+ Rl^JbCOii[2i)b] 

= I^kT ) ,,.0 ^ '"^ A ^ 

(Kl + Jl - Jl) - ?^^^:^Jb cos[2V^J cos[2V^,] = (4.54) 

This leaves only 9a with non-zero frequency. Therefore, critical points found in the 
[Nb, 0] polyads correspond to PO in the full phase space. 

Visualization of Critical Point PO In order to understand the critical points in 
a more intuitive manner, their corresponding PO are calculated in the Cartesian 
coordinate. This is done by assuming each normal mode to be a 2-dimensional har- 
monic oscillator with frequencies or uj^. Following the method in jHU], the normal 
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Cartesian coordinates are: 

X4 =[{Ka + K, + Ja- Jb)M]'/' C0S[(^, + + " ^b)/2] + 

[{Ka -Kb + Ja + cos[(^a - + + i^b)/2] (4.55a) 

[{K, -K, + Ja + Jfe)M]i/2 sin[(0, - ^6 + + ipb)/^] (4.55b) 

X5 =[(ira + - + Jb)/^^^]'^' COs[(^^„ + " + ^6)/2] + 

[(ir, -K,-J,- Jb)/uJ5?^' COs[{da -e,-^a- A) f^] (4.55c) 

t/5 =[(ira + i^b - + Jb)M]'/' Sin[(^, + eb-^a + A)/^]- 

[{Ka -Kh-Ja- Jb)I^^Y'^ M{Oa ' 9^ - i^a ' A)/2] (4.55d) 

(x4, ?/4, X5, 7/5) are related to the local Cartesian coordinates (xi, X2, 2/2) shown in 
Fig. ill by 

Xi = X4 + X5, 2/1=1/4 + 1/5 (4.56) 

X2 = -X4 + 2/2 = -2/4 + 2/5 (4.57) 




FIGURE 4.6 Cartesian bending coordinates, as defined in eqns. ()4.56|4.57p . 
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In order to visualize the resulting 3-dimensional bending motions, they are 
converted into animations with the modeling software Bryce 4 |l(J(Jj . The animations 
for Nh =22 are included on the accompanying CD-ROM in QuickTime 4 format |101j . 
Fig. 14.71 displays some of the still frames of these animations. The trajectories of L 
and Pre are both planar and resemble their namesakes in j79J. L has most of the 
amplitude of bending in one C-H bending oscillator (bender), and very little in the 
other. Pre has the two C-H benders at approximately equal amplitude but out of 
phase by |. For Orth, the two benders vibrate on planes orthogonal to each other, 
and reach their turning points in phase. Finally, in CR the two hydrogens rotate in 
ellipses (or circles when near Ja = 0) in opposite directions at the two ends of the 
C-C unit. 
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(d) (e) (fl 



FIGURE 4.7 Cartesian PO corresponding to critical points in [22, 0] polyad, produced 
from still frames in the computer animation. The carbon and hydrogen atoms are 
represented by the large and small spheres, respectively. The panels refer to (a) Trans; 
(b) Cis; (c) L; (d) Orth; (e) Pre; (f) CR critical points. 
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Immediately after each bifurcation, in Fig. 14.51 the nascent family of critical 
points remains close to Ja = ^Ka- The corresponding Cartesian trajectory resembles 
the "parent" normal mode motions. As the family migrates towards Ja/Ka =0, the 
Cartesian trajectory looks increasingly like the respective ideal motion in Fig. 14.71 
For example, the L trajectory initially resembles a slightly asymmetrical trans- bend. 
Then the imbalance of amplitude between the two C-H benders increases with A^";,. 
Finally, at =22 where Ja/Ka ~ 0, almost all the amplitude is in one of the two 
oscillators. 

The L and CR critical points are the global extrema points of the reduced 
Hamiltonian: L has the lowest energy while CR has the highest. This is apparent 
from an examination of the Vddi and Ve terms in ()4.28t I4.29|) reveals that both of 
them contribute most negatively to the Hamiltonian when [ipa, ipb) = (0, 0), and most 
positively when {ilJa,i^b) = (f , f)- 

In summary, critical points analysis of the [A*";,, 0] polyads yields four new 
families of critical points, in addition to the Trans and Cis normal families. The low- 
energy (EE) stability of Trans and Cis is destroyed in their first bifurcations. Due to 
the special condition Ki, = i = 0, the frequency of one cyclic angle 6b vanishes for 
all critical points. Therefore, the critical points correspond to PO in the full phase 
space. 
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4.3.3 Results of the [A^";,,^] Polyads 



Although the £ 7^ cases have not been explicitly considered by any other re- 
search groups, our analysis can be extended to arbitrary I values ^. Here we present 
critical points found in [A^{,,2], [A^^^, 6] and [A'^^, 10] polyad series as preliminary pre- 
dictions, and as a demonstration of this significant extension of our methods. These 
series correspond to = i/2 = 1, 3, 5, respectively. 

Similar to the i = case except using i^^end ^qn. ()4.19p as the Hamiltonian, 
two of the critical points equations are solved by holding (ipajipb) fixed at the four 
values in (j4.44p . However, it is no longer true the second of the two remaining 
equations 



has the trivial solution = as in the i = case. Instead, (|4.58( I4.59p must be 
solved simultaneously for the ( Ja, Jb) values. The solutions are numerically found by 
first transforming these equations into a polynomial form, then use the homotopy 
continuation package PHCpack |102j . This package has no limitation on the number 
and form of polynomials to solve, and can be used with no preliminary knowledge 
about the solutions. 

^Nevertheless, the validity of the effective Hamiltonian beyond £ = 2 is question- 
able due to a lack of supporting experimental data. 












(4.58) 



(4.59) 
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The JalKa and Jb/Ka values of the resulting critical points are plotted against 
Nb in Fig. 14.81 The solutions for each of the four [ipa-, iph) conditions are named in 
the same way as Table lOl In panels (a), (c) and (e), all families branch out of the 
Trans and Cis normal modes in a manner similar to the £ = results fFig. I4.5|l . From 
the three panels on the right, these new critical points are clearly not restricted to 
Jfo = 0. Instead, they diverge from J;, = with increasing £. 

The close resemblance between Fig. 14.51 and panel (a) of Fig. 14.81 is consistent 
with the observation that in DF spectra the £ = 0/2 states have nearly indistinguish- 
able intensity patterns. Were the dynamics vastly different for these two ^ series, 
the fractionation patterns of the bright states would be different and therefore dis- 
tinguishable. At still higher £, our analysis reveals interesting bifurcation structures, 
which awaits further interpretation. Even though the quantitative predictions here 
are limited by the validity of the effective Hamiltonian, we nevertheless expect the 
vibrational dynamics with sufficiently high I to become qualitatively different from 
the £ = case. 
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In spite of their similarity, the i = and i ^ critical points have different 

physical natures in the full phase space. In the former case they are periodic orbits 

and in the latter quasiperiodic orbits. The claim made by Jacobson et al. that "a 

fixed point in the abstract action-angle space . . . must lift to a periodic orbit in the 

physical coordinates of the molecule . . . that lies on a full dimension torus" |69j is 

accurate only when applied to i = 0. Generally, there exist as many independent 

frequencies as the number of additional constants of motion (polyad numbers). It is 

the critical points in reduced phase space, not the PO in the full phase space, that 

are the centers of phase space organization. 

4.4 Discussions 
4.4.1 Comparison with Other Studies 

Previously, the same Hbcnd Hamiltonian was analyzed by Jacobson et al. in 
a numerical study. They first visually examined sequences of SOS (of the 2 DOF 
reduced phase space) in order to identify stable fixed points at the center of regular 
regions. These fixed points are then followed while the energy is varied within the 
same polyad. A "family tree" of these fixed points is then built for [22, 0] [69j and 
[16, 0] |lU3j . In both cases, they found a local fixed point at the bottom and a counter 
rotator fixed point at the top energy end of the polyad. The fixed points were used 
to assign eigenstates based on the nodal pattern of the semiclassical wavef unctions. 

The local and counter rotator fixed points correspond to our critical points 
with the same names, because they occupy exactly the same place in the reduced 
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5 10 15 20 25 30 5 10 15 20 25 30 



FIGURE 4.8 Critical points and their bifurcations in [A*";,, I] polyads. The left panels 
display Ja/Ka vs. A^^, while the right panels display Jh/Ka vs. iVft values. Panels (a), 
(b): ^ = 2■ (c), (d): £ = 6; (e), (f): £=10. 
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phase space ^. There is another fixed point Jacobson et al. call the "M2 niode", 
which appears to be related to our Orth family of critical points with (EH) stability. 
Reproduced in Fig. 14.91 is the stable M2 PO trajectory they found in A^^ = 22. The 
motion qualitatively resembles panel (d) of our animation in Fig. 14.71 with the two 
C-H benders perpendicular but in phase with each other. Jacobson et al. observed 
numerically that this PO family becomes unstable below 14, 161 cm^^ within this 
polyad, while our Orth critical point at 14, 114 cm~^ has (EH) stability. Currently, 
we believe that the M2 family (both the stable and unstable segments) is formed by 
secondary PO surrounding the Orth critical point, which corresponds to a primary 

PO m- 



^Specifically, in [69] the local mode states are localized around {ipa, 4'b) = {0, 0), and 
the counter rotator mode states are around {n,7i). Their ipa and ipb are defined as 
twice ours. 
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FIGURE 4.9 The M2 trajectory in Fig. 3 of [^3, with A^^^, = 22 and E = 14, 561cm-^ 
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Independently, Champion et al. deduced the existence of L, CR and Orth 
modes as ideahzed cases using Lie algebraic analysis of the same bending Hamiltonian 
[73] . even though they did not perform a bifurcation analysis with variable polyad 
numbers. Moreover, eigenstates corresponding to the L and CR modes are located in 
ab initio calculations [Sni as well as numerical PO search on a molecular PES [7Hj . 

According to our analysis, the Pre family is bi-unstable (HH) and therefore 
expected to be surrounded by chaos. Thus, this family is not apparent in the study 
by Jacobson et al. , which relies on inspection of SOS for regular regions. Their 
SOS (reproduced below as Fig. 14.11)]) at nearby energies shows only a strongly chaotic 
region. 
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3ii/4 



FIGURE 4.10 SOS near the Pre critical point calculated by Jacobson et al. , and 
adapted from Fig. 1 of ^103. . The ipayi^b values are converted to the notation used in 
this thesis. Both sections have ipa fixed at ^. 
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4.4.2 Summary of Method 

Compared to the numerical procedure of Jacobson et ah , our analysis is more 
efficient in picking out major changes in the phase space structures of the 2 DOF 
reduced Hamiltonian here. In the former, a large number of SOS with different direc- 
tions and energy is examined in each polyad. In contrast, finding critical points in the 
reduced phase space involves only the solution of simultaneous analytic equations for 
a continuously varied polyad parameter. In the current case, almost all the significant 
observations from inspecting the SOS (specifically, those with manifestation in the 
quantum wavefunctions) can be obtained through four families of critical points. 

In the reduced Hamiltonian, time evolutions of the cyclic variable(s) are for- 
mally factored out and not considered as the essential part of dynamics. The critical 
points are determined in a uniform manner, regardless of the number (> 1) of polyad 
numbers. Hence, we are able to carry the same critical points analysis to both i = 
and i cases. At high i values, the latter cases display additional bifurcations that 
have not been observed previously. 

Moreover, the usage of SOS is almost exclusively reserved for 2 DOF systems. 
With 3 DOF, the phase space is 6-dimensional. A 2-dimensional surface does not have 
the correct dimensionality to divide it. Usage of 2-dimensional surfaces therefore 
is at best limited to short time scales, before the trajectories stop intersecting the 
dividing surface |lU4j . The generalized SOS is a 4-dimensional "hyper slice" whose 
visualization is not trivial JlU5j. Analytic detection (without visual inspection) of 
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dynamics becomes essential for analyzing systems with 3 or more DOF. The critical 
points analysis satisfied this requirement. The number of equations scales linearly 
with the DOF. In addition, all critical points are found in a comprehensive manner 
regardless of their stability property. The resulting critical points then enable a guided 
exploration of the phase space using other techniques. 

4.5 Effect of Single DD-I and Single I Resonances 

In this section, we propose an interpretation of the bifurcations in the full 
4 DOF iffcend in terms of its two strong resonances ( namely the DD-I and the ^ 
resonances) taken individually. With either one of these resonances, the [Nb, £] polyad 
is further separated into subpolyads. The DD-I resonance subpolyads are ahgned 
along the columns in Fig. \A.'S\ while the £ resonance subpolyads are along the rows. 

The corresponding single-resonance Hamiltonians 

Hddi = Ho + Vddi (4.60) 
H, = Ho + V, (4.61) 

can be analyzed with the method described in § 3.1 with m : n = 2 : 2. The Huui 
Hamiltonian is parameterized by two reduced phase space variables {Ja.ipa) according 
to eqns. (j4.18|) . Within each DD-I subpolyad, these parameters define a DD-I PPS. 
Similarly, the Hi Hamiltonian is reduced to ( J^, ipb) space, which defines an £ PPS 
for each £ subpolyad. For simplicity, we only consider the [A*";,, 0] polyads without the 
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high-order coefficients yijkjr^^^ and r^^^. Then each DD-I PPS has radius Ka — \Jb\ 
while each i PPS has radius Ka — \Ja\- 

Fig. 14.11] presents the zero-order DD-I and £ PPS, and the semiclassical tra- 
jectories on them. With only the Hq term, the Hamiltonian does not depend on ifja or 
ipb. All semiclassical trajectories of Hq can be regarded as the "composition" between 
one DD-I trajectory and one i trajectory, as labeled in the figure. The critical points 
on the PPS are where 

dHnnj dHnnj dHp dHf 

In the zero-order case, the only critical points are the north and south poles of each 
PPS. 

The DD-I and i PPS for polyads [4,0], [12,0] and [20,0] are presented in 
Figs. 14.12114.1^ Note that the two sets of spheres do not coexist within the same 
Hamiltonian. The zero-order quantum numbers on these figures are no longer mean- 
ingful in labeling trajectories in Hddi and Hi, since n4,n5 (or ^4,^5) are mixed by 
the resonance term. As the DD-I and £ PPS do not qualitatively change across the 
subpolyads with same Nb, in the latter two figures the number of PPS is reduced to 
make the figures legible. 
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The HoDi System This system had been considered previously by Rose and Kell- 
man with shghtly different parameters but the same quahtative results |79j. A com- 
parison across Figs. 14.1^4.141 shows that as Nh increases, two new families of critical 
points emerge in two bifurcations on the DD-I PPS. Between Ni, =6-8 (Figs. 14.121 
I4.13|) . two local mode critical points L^iji are born at the north pole of the PPS onto 
the great circle {ipa = 0,7r). Then between Nf, =14-6 fFigs. I4.14|) two preces- 
sional mode critical points Pddi emerge at the south pole onto another great circle 
ii^a = f , Both Lddi and Pddi migrate towards the equator of the PPS (where 
Ja =0) with increasing A^";,. 

The migration of the Lddi and Pddi is depicted in Fig. 14.151 which plots these 
critical points in the same manner as Fig. 14.51 Similar to Fig. 14. 5| here both critical 
points in Hddi are born at Ja/Ka = ±1 and migrate towards the equator, finally 
reaching Ja/ Ka = at about Nb = 20. 



117 




FIGURE 4.15 Bifurcation of critical points in Hddi for the central column subpolyads 

(4 = 4 = 0). 
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Fig. 14.161 summarizes the structure of the central DD-I PPS {Nh =4-20) on 
the 2 : 2 catastrophe map. The control parameters Sddi, P'ddi defined in eqns. 
(4.1-4.2) of pHl- P' characterizes the strength of the 2:2 resonance, while 6 refiects 
the detuning between the two coupled frequencies. Here, the resonance strength S45 
is a constant. The trans- and cis- bending frequencies, starting as 004^ < cus, are tuned 
towards each other as X55 < < X44. The representative points on Fig. 14.161 cross two 
zone boundaries as Nf, increases: first the trans- normal critical point (north pole) 
is destabilized in a bifurcation, then the cis- normal critical point is destabilized by 
another bifurcation at higher A*";,. 

Ref. jini assumes a planar model of C2H2 bending (see Fig. 1 thereof), and 
the Lddi and Pddi critical points have the same Cartesian motion as the L and Pre 
in Fig. 14.71 These two modes can be interpreted simply as follows. The angle i/ja 
is defined as the relative phase between the trans- and cis- normal mode oscillators. 
When exactly in resonance, these two oscillators have the same frequency. As shown 
in Fig. 14.181 (a) and (b), superimposing these two oscillators on the same plane with 
relative phases or vr results in only one C-H bender being excited (L). Changing the 
relative phase to | or ^ causes the two C-H benders to be out of phase by | or ^ 
(Pre). 
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FIGURE 4.16 Catastrophe map of Hddi- The representative points are for the central 
DD-I PPS with Ni, = 4-20. N, L and P denote the normal, local and precessional 
mode critical points. The unstable critical points of a given region are underlined. 
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FIGURE 4.17 Catastrophe map of H^. The zone structures are the same as Fig. 14.161 
At all values, the representative points of the central ^ PPS are located at the 
same point. 
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The System Unlike Hddi, the Hamiltonian Hi satisfies a special degeneracy 
condition. The quantum ZOS \n^^,nl^) and \n^^'^ ,n^^^) always have the same energy. 
This is a consequence of the time reversal symmetry, when the directions of both 
angular momenta £4, ^5 are reversed |1U6 ]. In other words, the two "frequencies" 
coupled by Vi are always in exact resonance. Any finite Vg term will induce a non- 
local bifurcation of critical points on the zero-order system jin]. Such a bifurcation 
is indeed shown when we compare the i PPS between Figs. 14.111 and 14.121 When 
Ve is turned on, both poles of the i sphere become unstable. The following two new 
families are stable critical points: 

Lf. Jb = 0, z^b = 0,7r (4.63) 

P^: Jb = 0, V'6=2'Y (4-64) 

Because they are both 2 : 2 type resonances, Hijijj and Hi have the same zone 
structure on their catastrophe maps. The latter is shown in Fig. 14.171 The coupling 
strength is taken to be r^^ (which is a constant), and f3' = 0. On the catastrophe 
map, the representative points for the central i subpolyads in Figs. I4.l21l4.14l share a 
single location for all values. 

Formally, the Li and Pe modes are like the local and precessional modes in 
Hddi, although the angles ipa and ipb have very different meanings. While ipa is the 
relative phase angle between planar trans- and cis- motions, the physical interpreta- 
tion of ipb is rather complex. When £4 = £5 = 0, ipb corresponds to the dihedral angle 
between the two planar normal C-H benders. 
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FIGURE 4.18 New bending modes as superposit ions of the normal trans- and cis- 
modes. The coordinates are the same as in Fig. 14.61 For quahtative purposes, the 
trans- and cis- vibrations are assumed to have the same frequency and amphtude. 
In panels (a-d), the relative (phase angle, dihedral angle) between the normal mode 
oscillators are in the same order as in eqn. ()4.44|) . Superimposing trans- (first column) 
and cis- (second column) motion result in the four new types of motion found through 
critical points analysis of the full Hamiltonian Hbend- 



123 

Connection with Critical Points in Hf,end A comparison between the critical points 
in H£)£)j and Hi and those found in the full Hbend of § 4.4 reveals the following 
connections: 

1. The 16 conditions of (ipayipb) for critical points in Hbend (Table can be 
produced from combining the great circles on DD-I and i PPS. Table can 
thus be regarded as a set of generalized great circles in H^end- 

2. DD-I resonance induces bifurcation of critical points in ( J^, ipa) coordinates with 
increasing Ni^. The new \^ddi and Pdd/ then migrate towards = 0. The £- 
resonance, on the other hand, induces and P^ critical points, both located 
at Jf, = 0, as soon as Vi is included. There is no additional bifurcation in Hi as 
iVft increases. All these observations are consistent with the results in Fig. 14.51 

3. The four types of non-normal PO in Fig. 14. 71 can be produced by superimposing 
equal amounts of trans- and cis- vibrations at the same frequency, with rela- 
tive phase angle ipa and dihedral angle iph- This superposition is illustrated in 
Fig. 14.181 adding the trans- and cis- bends with (relative phase, dihedral angle) 
at (0,0) and (0, ^) result in the L and Orth type vibrations (panels a and b), 
while (f ,0) and (|, |) result in the Pre and CR type vibrations (panels c and 
d). 

4. Finally, the critical points in Hbend can be regarded as formed from one critical 
point in Hodi and one in Hi, in the manner listed in Table WM The bifurca- 
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tions along Nf, are caused by DD-I, while the nature of the four new modes is 
determined by both DD-I and i resonances. 



TABLE 4.4 Proposed composition of critical points in H^endjiom those in H^uj and 
Hi. For the notations in columns 1 and 3, refer to Figs. 14. 12114. lil 



Hddi and He 


Hhend 


Hddi and He 


Hbend 


^DDi, Le 


Trans 


hoDi, Le 


L 


Ndd/, Pe 


Trans 


Pddi, Le 


Orth 


Sddi, Le 


Cis 


Ldd/, Pe 


Pre 


S_DZ)/5 Pi 


Cis 


Pdd/, Pe 


CR 



This composition is only a qualitative one. Most notably, it fails in predicting the 
stability of critical points in Hbend- Although all new critical points in Hddi and 
He are stable, those in Hbend exhibit three different types of linear stability. Only 
two families, the L and CR, are bi-stable (EE). Another discrepancy is that while 
Hbend has 4 distinct bifurcations, simply considering the single resonance systems 
according to Fig. 14.151 would suggest that L and Pre, Orth and CR are born in only 
two bifurcations. These differences are likely to be caused by the fact that in the 
full Hamiltonian, the two directions ( Ja, ^a) and {Jb, ipb) are strongly coupled to each 
other, instead of forming independent subsystems Hddi + He- 

4.6 Quantum Survival Probabilities 

Due to its resemblance to the transition state of Fig. 14.11 the local bending 
mode is expected to play an important role in the acetylene-vinylidene isomerization 
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dynamics. Recently Carter et al. performed a Car-Parrinello type calculation on this 
system [107j. Dozens of classical trajectories are integrated, with the forces at each 
step obtained from ab initio calculations. Surprisingly, many trajectories go back and 
forth across the barrier many times, before settling in the acetylene potential well. 
The authors give a kinematic explanation: The time interval between the hydrogen 
crossing the barrier and swinging back is not long enough for the energy in the 
reaction coordinate to dissipate effectively. This interpretation reconciles an existing 
discrepancy of vinylidene lifetime - 0.04-4.6 ps in |108|109j and 3.5 us in |110j . The 
picosecond timescale is believed to be that of the initial decay of vinylidene, while 
the microsecond one is the vinylidene lifetime averaged over many recrossings. 

Such a "recurring state" is necessarily decoupled from the rest of the vibra- 
tional manifold. This is supported by the observation of Levin et al. In their 
Coulomb explosion experiment, while the vinylidene molecules have energy well above 
the reaction barrier, the estimated dilution factor is only ^ 0.5, indicating that this is 
coupled to about one other state. In another independent study, Schork and Koppel 
compared the intrinsic lifetime of vinylidene to the local density of acetylene vibra- 
tional states, and concluded that extensive IVR is unlikely, at least for the lowest 
vibrational state of vinylidene Srivastava et al. also suggested the acetylene- 
vinylidene isomerization is going to deviate significantly from the RRKM limit, due 
to the relatively low density of states |112j . 

The survival probability P(t) has been an important tool in characterizing 
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the dynamics of quantum states |74|113j . For a quantum state written as an 
expansion in the eigenstate basis 

\^it)) = J2c^e-'''^'/^\<P^) (4.65) 
P{t) is defined as the overlap between |^I/(0)) and |\I/(t)): 

P(t) =|(M/(t)|M/(0))r = I (^X:(0.|c:e^^»*/^j (y1c,\^j) 
^\c,\^cos[Eit/h]^ + (^\cj\^sm[Ejt/h] 

= ^\ci\^ + '^ \ci\'^\cj\^cos[{E^- Ej)t/h] (4.66) 

The survival probability of an eigenstate is trivial: P(t) =1. For other initial 
states, the initial decay from P(0) =1 is dominated by how many states are directly 
coupled to 1^) jll4j . Oscillations (quantum beats) in P{t) at intermediate time 
scales describe the usually partial recurrences of the initial state. Finally, the long 
time average ^ also known as the dilution factor jS], gives an estimate of the 
number of states participating in the IVR of the initial state. 

P{t) is the quantum analogue of the classical autocorrelation function |115j . 
When it remains near unity for a sufficiently long time, then is strongly localized 
in some representation (as opposed to spreading over all space). In the statistical 
(RRKM) limit, P{t) quickly decays to a value close to the local density of states. 

In § 3.3.3, we hypothesized that semiclassical wavefunctions are likely to lo- 
calize near the critical points that are also extremum points in the reduced classical 
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phase space. To test this claim, P{t) of states corresponding to Trans, Cis, L and 
CR overtone states are calculated for polyads A*";, = 2 to 20. While the Trans and 
Cis states are the normal ZOS |n4,0°) and |0'',?t,5), the L and CR states have to be 
constructed. Here we use the method described by Field et al. in f7^ (which also 
contains calculation of P{t) for selected L states). The L and CR states are defined 

7 

as : 

|Ln) = ^(«L + «L)(«l + 4,)^|0°,0°) (4.67) 
|CRn) = ^(aL + aL)(«L-«L)^|0°,0°) (4.68) 

These are the "perfect" L and CR states in the sense of containing equal amounts of 
trans- and cis- components. 

Fig. 14.191 shows the P{t) values for the first 4 picoseconds. Similar results had 
been obtained by Jacobson et al. in Fig. 3 in (Sn]. In panels (a,b), up to Njy =6, 
the Trans and Cis states have P{t) ^ 1, indicating that these overtone states are 
well localized near the respective critical points. At A^^ =10, the periodic oscillations 
become stronger, but the recurrences still reach close to 1. The stability appears to 
be lost for the Trans state at A";, =10 and for the Cis state at A^^ =14. These changes 
occur soon after the first bifurcations of the Trans and Cis critical points in Fig. 14. 5| at 
Nfy =8 and 10 respectively. However, even at A";, =22 their dilution factors are still 0.4 
and 0.3 for 122°, 0°) and |0°,22°), respectively. These two states are now coupled to 

derivation for the L state can also be found in § 5.3.2. of |116j . 



128 

71 other ZOS states in the polyad. Had the IVR been purely statistical, the dilution 
factor would have been an order of magnitude smaller than the actual dilution 
factors. This suggests the vibrational dynamics are far from the RRKM limit, even 
after then normal modes have been destabilized. 

In panels (c) and (d) of Fig. I4.19[ before A'';, = 14 the P(t) evolution of L and 
CR states remain strongly oscillating. The almost sinusoidal oscillation between 
and 1 at A*";, =2 is due to the fact that each of these states is coupled to another L 
or CR states that is degenerate to it. Between Nb =18-22, both L and CR states 
have P{t) oscillating slightly under unity. On the bifurcation diagram of Fig. 14.51 
this corresponds to (J^ = 0) which is where the (EE)-type critical points L and CR 
approach their "prefect" shape in the Cartesian coordinates. Hence, quantum wave- 
functions localize around these now stable modes. Similar conclusions were reached 
by Jacobson et al. through a visual match between the classical PO and semiclassical 
wavefunctions |69|103j . 
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In conclusion, the states in Fig. 14.191 with near unity P{t) are in agreement 
with the (EE)-type critical points found in the classical analysis. P(t) of L states 
in Nb =22 shows a lack of significant IVR, which is similar to the behavior of the 
isomerization states discussed in |llUj . The only difference is that the pure bending 
acetylene states presumably are still below the isomerization barrier. With additional 
vibrational energy in the stretching DOF, these states or the superposition of several 
such states will be capable of going back and forth across the barrier many times. 

4.7 Summary and Conclusion 

The critical points analysis is performed on the C2H2 pure bending effective 
Hamiltonian. In the [A*";,, 0] polyads, 4 new families of critical points (L, Orth, Pre and 
CR) are born out of the normal Trans and Cis critical points in distinct bifurcations as 
Nb is increased. The bifurcation points where the new families are born correspond 
to qualitative changes in the classical phase space structure. Computer-generated 
animations give visual insight into the nature of their motions in Cartesian space. 
Similar bifurcation structure is obtained for the [A'';,, i] polyads with £ = 2,6 and 10. 

Three of these new families are consistent with the results of other researchers, 
who had used more elaborate methods. The bi-unstable Pre family, on the other 
hand, can only be uncovered through an explicit search of critical points like ours, as 
presented in this thesis. 

Separate consideration of the DD-I or i resonances qualitatively accounts for 
the origin and nature of these four new modes. Using the method of § 3.1, the 
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dynamics induced by DD-I or I resonance alone is analyzed. Combining the critical 
points in the single-resonance Hamiltonians yields the same types of motion as those 
obtained as critical points in the full bending Hamiltonian. 

The calculated quantum survival probability shows the break down of the 
normal modes description at intermediate polyad, as well as the emergence of L and 
CR modes as new stable modes of vibration at Ni^ =22. These results demonstrate 
that the classical phase space structure is indeed reflected in the dynamics of the 
corresponding quantum system. 
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CHAPTER V 

BIFURCATION ANALYSIS OF C2H2 STRETCH-BEND 
5.1 Introduction 

Beyond the pure bending subsystem, the next step is extending the critical 
point analysis to all 7 vibrational DOF (4 bending and 3 stretching) of C2H2 explicitly. 
In such a multidimensional system, an analytic detection method would be superior 
to both visual inspection and numerical search. 

5.1.1 The Effective Hamiltonian 

Three stretch-bend effective Hamiltonians of C2H2 have been published up to 
date |fj:jlll7lll8j . 

In this chapter, we use the first Hamiltonian of Herman et al. , because none 
of the other Hamiltonians has been studied theoretically |8U|82j . The Hamiltonian 
in 1221 contains 8 resonances: one Darling-Dennison type -ft'n/ss resonance coupling 
between the two normal C-H stretch modes, the -ft'44/55 (known as DD-I in Chapter 
4) between the two normal bending modes, the £ resonance r\^, and -ft'3/245, 
-^1/255, -^14/355 -^33/1244 resouauccs which couple between the stretch and bend DOF. 
This coupling structure is illustrated in Fig. 15.11 In a typical stretch-bend polyad, 
the normal ZOS |ni, ^2, n3, n4*, ng^) are all coupled by a complex web of resonances 
except for two special cases. The pure bending polyads analyzed in Chapter 4 form 
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an isolated subsystem with polyad number Ng — 0. The C-C stretch overtones 
|0, n2, 0, O'', 0°) are not coupled by any resonance; therefore they are eigenstates of the 
Hamiltonian. 
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Vi Sym. C-H 
Stretch 

V3 Asym. C-H 
Stretch 

V2 C-C Stretch 
V4 Trans-Bend 
V5 Cis-Bend 



4 Vibrational 
Angular Momenta 



Stretch DD 



Bend DD (K44/55) 



Stretch-Bend 



K3/245 , Ki/244 , 
Ki/255 , K14/35, 
K33/1244 



C resonance (r45) 



FIGURE 5.1 Resonance couplings among the stretch-bend normal modes, according 
to the Hamiltonian of ^92jj. On the left are the normal modes of Fig. 14.21 



135 

In producing this Hamiltonians (as well as the other later ones), the fitting 
is based on only a small part of the states predicted by the polyad model. A high- 
lying polyad typically contains dozens of states, out of which only a few have been 
experimentally detected and included in the fit. This allows for some uncertainty in 
the high-order coefficients in the fit. Especially lacking are states that contain both 
stretching (A^^ > 1) excitation and high bending excitation (n4 + ^5 > 12). 

Using Heisenberg's Correspondence Principle in eqn. (|2.11|) . the 7 DOF quan- 
tum Hamiltonian in [^21 is transformed to a classical one with 14 action-angle vari- 
ables: (tj, (pi) for i = 1 — 5 and {kj, xj) for j = 4 — 5. The actions are related to the 
zero-order quantum numbers by 

ri = ni + ^ fori = 1,2, 3 (5.1) 

Ti = ni + 1 for i = 4, 5 (5.2) 

Kj = ij for j = 4, 5 (5.3) 

To simplify our analysis, the last five high-order parameters in Table 15.11 
?/244, -^"33/1244, ki, r445, r545, are ignored in the rest of this Chapter. The classical Hamil- 
tonian IHOl then has the following form: 

Hsb =Ho + Ha (5.4) 

with 

55 5 

Ho = '^UJiTi+ ^ XijTiTj + ^ gijKiKj (5.5) 
«=1 'i;j=hi<j i,j=i;i<j 
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rir3Cos[2((/)i - (^3)] 



+ 



K- 



3/245 



V'^2T3{\/ (t4 - K4)(7"5 + K5) COs[(/)2 - (/)3 + (^4 " X4 + (/"S + Xs] 



+ 



+ V (r4 + K4)(t5 - K5) COs[</>2 - 03 + '^4 + X4 + '/'5 - Xs]} 



244 



\ItiT2{tI - k\) cos [01 - 02 - 204] 



+ 



1/255 



nT2(T-| - K5) COS [01 



205 



K 



+ 



(r4 - «;4)(7-5 - K5) cos [01 - 03 + 04 - X4 - 05 + X5] 

5 -X5]} 



+ V (r4 + K4)(r5 + K5) COS[01 - 03 + 04 + X4 



+ 



+ 



?'45 



pV(^|-'^I)(^|-'^i)cOs[2(04 - 05)] 
\/ (t-4 - /^4)(^5 - '^i) C0S[2(X4 - X5)] 



(5.6) 



The values of the parameters in Hgi, are hsted in Table 15.11 ^ The 3 polyad 
numbers correspond to 3 classical constants of motion: 

15 

P =5ri + 3r2 + ^t^ + + = Nt + — (5.7) 

R=T^+T2 + T^ = N, + ^ (5.8) 
L =K4 + K5=e (5.9) 

5.1.2 Overview of Existing Studies 

To our knowledge, the stretch-bend effective Hamiltonians has never been 
analyzed with all the resonances. We are aware of only three relevant studies: 

^It was found later that the -ft'14/35 value appears as 29.044 in [^21 and as 29.944 
in jHO]- This appears as a minor discrepancy compared to its later revision to 15.66 
in dnj and then to 16.614 in [TTR] . 
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TABLE 5.1 C2H2 stretch-bend effective Hamiltonian from [22] and published in [HHj 
for the classical Hamiltonian. The parameters are in units of cm~^. 





3501.537 


2^15 


-10.09 


2^45 


-2.311 


-^14/35 


29.944 


Ll)2 


2013.425 


2^22 


-7.802 


2^55 


-2.492 


-^44/55 


-12.909 


CJ3 


3417.644 


2^23 


-5.882 


5'44 


0.4181 


^45 


-6.09 




621.692 


2^24 


-12.841 


9ib 


6.603 








746.773 


2^25 


-1.829 


955 


3.676 


1/244 


0.1522 


Xii 


-24.758 


2^33 


-27.483 


Kn/33 


-102.816 


-^33/1244 


6.38 


X12 


-11.199 


2^34 


-10.617 


-^3/245 


-16.698 




-1.315 




-103.386 


2^35 


-8.676 


-^1/244 


6.379 


^AA5 


0.1255 


Xm 


-12.98 


X44 


3.595 


-^1/255 


6.379 


'f' 455 


-0.225 



• Pals and Gaspard ^SOJ investigated the recurrences of the classical trajectories. 
Although the study did include the stretching DOF, the authors were mostly 
focused on the pure bending subsystem. 

• Hasegawa and Someda jH21 analyzed the quantum dynamics using a perturba- 
tive method. The focus was the short-time evolution of 3 types of quantum 
ZOS: |0,0,0,nO,0°), |1, 3, 0, 6°, 0°) and |0, 0, ng, 0°, 0°). 

• Kellman et al. investigated the planar system (£4 = ^5 = 0) with 3 resonances: 
i^ii/33, _K'3/245 and -ft'44/55 in a diabatic correlation approach |119j . A series of 
states termed the "primary subpolyad" are identified, which carries most of 
the intensity in the experimental spectra. The states are then fit to a single- 
resonance Hamiltonian in an effort to interpret the spectral patterns. 
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Due to the dimensionality of the problem, none of these studies explicitly 
considered all the major resonances in the effective Hamiltonian. Since the approach 
presented in this thesis is designed to be dimensionality-independent, we believe it is 
more suitable for analyzing the stretch-bend acetylene system. 

5.2 Preliminary Considerations 

In Chapter 4 it was shown that the normal bending modes are destabilized 
at increased internal energy (polyad numbers) by the resonance couplings, and new 
stable bending modes are born in the bifurcations. In the stretch-bend system, the 
PO search by Prosmiti and Farantos also indicates that the highly excited stretch- 
bend system retains some regularity |12Uj . Therefore one can expect a similar to the 
pure bending case to exist in the stretch-bend dynamics, provided there are sufficient 
excitation and couphng. 

5.2.1 The Stretch Overtone Polyads 

The critical points analysis starts with solving for all critical points in the 4 
DOF (7 DOF - 3 polyad numbers) reduced classical Hamiltonian. In obtaining the 
preliminary results, we focus on a single series of polyads of interest, as opposed to 
varying all three polyad numbers independently. 

Of the five normal modes of C2ll2 depicted in Fig. 14.21 the pure bending 
subsystem was analyzed in Chapter 4, and the C-C stretch {T2) overtones form an 
isolated subsystem. This leaves the obvious question: what could happen to the 
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two C-H normal stretches, as they are excited to higher energy? Although the C- 
H stretching dynamics has been actively studied using two-mode models |7n|, their 
coupling to the other vibrational modes is complex and poorly understood. This 
motivates us to investigate the following polyad series, which includes the symmetric 
and antisymmetric C-H stretch overtones ^: 

{ri, r2, ra, r4, n, k^, k^} = {n, 0, 0, 0, 0, 0, 0}, {0, 0, n, 0, 0, 0, 0} (5.10) 

The polyads containing them are found by substituting eqn. (j5.1U|) into eqns. (|5.7I - 
15. 9|) to give 

P = 5R, L = i = (5.11) 

These [P, R, i] = [5R, R, 0] polyads are henceforward referred to as the stretch over- 
tone polyads. The upper limit of R in this study is set at 8, as the effective Hamil- 
tonian in |117j includes up to 6 quanta of C-H overtones excitation {R = 7.5). 
5.2.2 Stability of the Normal C-H Stretch Overtones 

The first consideration is the stability of the normal C-H stretch overtone 
states. For this purpose, the classical trajectories very close to the normal mode 
overtones in eqn. ()5.10|) are integrated. The overtone is classically stable if the devi- 
ation of nearby trajectories remains small. The onset of large-amplitude oscillations 
indicates the destabilization of the corresponding overtone trajectory. 

^Here in the purely classical analysis, the quantization requirement that the action 
in each mode exceed the zero-order energy (see eqns. 15.11 15. 2j) will be ignored. 
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The results are displayed in Fig. 15. 2[ In panel (a), the symmetric stretch 
overtone ti remains stable to at least ti = 6.5 {rii = 6). In contrast, the antisymmetric 
Ts overtone becomes unstable at as low as rs = 1.5 (na = 1). Two resonances may be 
responsible for the destabilization of t^, overtones: -ft'n/as and -^'3/245- 
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FIGURE 5.2 Classical stability of C-H normal stretch overtones. Panel (a) and (b) 
are the ti and T3 values of trajectories very close to the respective overtone conditions 
in eqn. 
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Single -ft'ii/as Resonance The single i^n/sa resonance system 

^^11/33 = ^^0 + ^^nrg cos[2(0i - 03)] (5.12) 

is analyzed for the stretch overtone polyads using the formulation of § 3.1, with 
m : n = 2 : 2. After the other actions T2,T4,t^, n^, are set to be zero, (Ti,r3) 
form a two-mode system. The 2 DOF reduced phase space is described be canonical 
variables I^, V'- 

/. = ^Y^, ^ = 01-03 (5.13) 

The constant energy contours (at arbitrary energy) of i^n/ss in the {Iz,ip) 
space are plotted in Fig. 15.31 Each of these panels is the Mecartor projection of the 
corresponding PPS, as in panel (a) of Fig. 13. 11 The top and bottom of Iz represent the 
Ti and overtones, respectively. The contours display the same qualitative features 
as semiclassical trajectories of the eigenstates. The critical points (dark dots in panel 
a) at = 0, TT and 27r are the stable (E) local C-H stretch critical points ^. The 
local model of C-H stretch has been previously used in studying the two coupled C-H 
stretch oscillators |121ll22j . 

The curves at the top of all panels of Fig. l5.3l that run across [0, 2tt] indicate the 
Ti overtone is a stable critical point. The overtone, in contrast, becomes unstable 
on the PPS from ts = 2 onward. These qualitative observations verify the stabilities 
shown in Fig. 15.21 

^The points at and 2tt are identified with each other. They are also energetically 
degenerate with the one at vr, as the two local C-H stretches are classically equivalent. 



143 



(a)R=2 (h)n=4 




FIGURE 5.3 Reduced phase space of -f/^11/33 Hamiltonian, and the constant energy 
contours. In panels (a)-(d) i? = 2,4,6 and 8 respectively. The dashed line in panel 
(a) labels the separatrix crossing the unstable Ts overtones (the = — 1 line). The 
black dots at ip = 0,7i and 2n indicate the stable local mode critical point. 
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Single -ft'3/245 Resonance The single -ft'3/245 resonance system 

7^3/245 I 

^^3/245 = ^^0 H ^ V (-r4 - K4)(r5 + Ats) COs[(/)2 - (/-S + (/>4 " X4 + (/>5 + Xs] 

+ -\/(r4 + K4)(r5 - K5) cos[(/>2 - 03 + </'4 + X4 + (/"s - X5]} (5.14) 

was formally analyzed by Rose and Kellman in [S3] ■ However, they failed in perform- 
ing a calculation using actual parameters of C2H2. Here we perform the calculation 
for the stretch overtone polyads, with actions ri = /t4 = K5 = ^. The r3 overtones 
are then contained in the action space with 

T2 = ta = T5 = R- T3 (5.15) 

Because the derivation in failed to use the proper di for the doubly-degenerate 
bends, the canonical variables (J, (p) ^ of the 1 DOF reduced phase space are redefined 
as: 

J = 14X2 + T3- 7X4 - 7X5, 95 = 03 - 02 - 04 - 05 (5.16) 

Fig. 15.41 displavs the arbitrary constant energy contours in {J,(p) space. The 
top of each panel (maximum J) corresponds to the X3 overtones and the bottom to 

X2 = X4 = X5 = i?, n = 

(5.17) 

^In jHS], the elimination of from consideration is not rigorous, since the -R'3/245 
resonance does couple £4,^5, and therefore ^4,^5. 

^The notations are intentionally different from those used in ifii/33 above to avoid 
confusion. 
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In panel (a), the resonance induces a bifurcation of the T3 overtone. In this 

process, a stable critical point is created ai ip = ir. However, the overtone itself 
remains stable at higher R (panels b-d). 




FIGURE 5.4 Reduced phase space of ^^3/245 Hamiltonian, and the constant energy 
contours. In panels (a)-(d) R— 2, A, 6 and 8. 
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From the above separate consideration of -ffii/33 and -^3/245, the overtones 

are first destabihzed by the -ft'11/33 resonance below R = 2, and the local C-H stretch 

becomes the stable mode. 

5.2.3 Effect of Stretch-Bend Resonances 

After its creation, the local C-H stretch may be further perturbed by the 
stretch-bend resonances to couple with the bending DOF. Of the four stretch-bend 
resonances (-ft'3/245, -f^i/255 and Ki^/^^), it is unlikely that they are equally 

important at all parts of a stretch overtone polyad. We use a simplified version 
of Chirikov analysis |123j in order to estimate of which resonances are important 
in affecting the local C-H stretch. When K4, K5 are left out of consideration, the 
"strength" of each of these 4 resonances as well as i^ii/33 can be measured by how 
far the zero-order frequencies Ui = dHQ/dri are tuned towards the integer ratio rrii 
corresponding to an exact resonance condition. At the exact integer ratio, Ui satisfy 

PI: 

i 

The vector with integer components {mi, m2, ms, 7714, ms} corresponds to the reso- 
nance vectors mentioned in § 2.1. In this case they are: 

Kn/33 :{2, 0, -2, 0, 0}, 7^3/245 :{0, -1, 1, -1, -1} 

i^i/244 :{1, -1, 0, -2, 0}, i^i/255 :{1, -1, 0, 0, -2} 

7^14/35 :{1, 0,-1, 1,-1} 
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For each resonance, eqn. ()5.18|) defines a 4-dimensional hypersurface in the 
T2, T3, T4, Ts} space. Because of the conservation of P and R, only 3 of the actions 
can be independently varied within a polyad. It is therefore possible to represent these 
hypersurfaces as 2-dimensional surfaces, which are referred to as resonance planes in 
a 3-dimensional volume (Ji, J2, J3) 

Ji=n-T, e[-R,R] (5.19) 

J2=T-2 e[0,R] (5.20) 

J3=n-n e[-2R,2R\ (5.21) 

This volume is further constrained by the requirement that all the r, should be non- 
negative. 

Fig. 15.51 depicts the resonance planes in {Ji, J2, -h)- In panel (a), at i? = 4 
the C-H stretch system (thick line) is separately in contact with -ft'n/ss (red) and 

(magenta) resonance planes at two places. The -^^1/244 resonance only perturbs 
but does not destabilize the ri overtones, according to Fig. 15.21 In panels (b) and 
(c), the C-H stretch system first interacts with the -^^11/33 resonance plane, which is 
next intersected by the -ft'1/244 plane. In all 3 panels the -ft'3/245 and A'1/255 planes are 
located at another side of the (Ji, J2, J3) space, and can only interact with the C-H 
stretch system via i^n/ss and other resonance planes. Because the -ft'14/35 resonance 
term vanishes in the absence of bending excitation, the -ft'14/35 its plane (navy) also 
cannot directly interact with the C-H stretch system. 
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5.3 Critical Points Analysis 



Based on the results of § 5.2.2-5.2.3, in this subsection we compute critical 
points in the following three cases: 

1. The (tijTs) subsystem with fTn/ss resonance; 

2. The (ti, r2, Ts, subsystem with both Knj^^, and ^1/244 resonances; 

3. The full Hsb of eqn. ()5.6|) with all 7 resonance couplings. 



Using the procedure described in Appendix A, a canonical transformation is se- 
lected that expresses iJ^b in the following new action-angle variables: 3 pairs of trivial 
action-angle variables (P, i?, L, Op, 9r, 9l) for the polyad numbers (defined in I5.7ll3^ 
and their conjugate angles, and 4 pairs of non-trivial ones (/i, I2, 13, h, "^i, "^2, ^3, ^4) 
which span the 8-dimensional reduced phase space. The angles \l'i — \&4 are chosen to 
correspond to the Kii/2,3, -ft'1/244, -^44/55 and i resonances, respectively. The simplest 
canonical transformation we have found so far is: 



5.3.1 Computational Details 



h =12ri + 8x2 + llrs + 2u + 2n + 2^4 + 2^5, 



^'l =01 - 03 



(5.22) 



I2 =6ri + 3r+6r3 + n + + k,4 + k,^, 



^2 =01 - 02 - 204 (5.23) 



-^3 =6ri + 4x2 + 6x3 + T4 + K4 + K5, 



*3 =04 - 05 



(5.24) 



h = - 



12ri — 8r2 — 12x3 ~ 2r4 — 2t^ — — 



^4 =X4 - X5 



(5.25) 
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FIGURE 5.5 Frequency resonance planes in the stretch overtone polyads. Panels (a), 
(b) and (c) display the cases with i? = 4, 6 and 8, respectively. The color coding 
is: Kii/33 (red), i^3/245 (turquoise), i^i/244 (magenta), Ki/255 (green), i^M/as (navy). 
The C-H stretch system with T2 = T4 = = is indicated by the thick black line in 
each panel. 
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is the relative phase angle between the symmetric ti and antisymmetric oscil- 
lators. \E'3 and \E'4 are identical to ipa and ipb in Chapter 4, respectively. 
Critical points in the reduced phase space are defined by: 

(5.26) 
(5.27) 

In eqns. (|5.26|) . only appear in the form of cos functions. Similarly to the 

treatment of the pure bending subsystem (§ 4.3.1), a sufficient condition for them to 

be simultaneously satisfied is: 

sin[2\E'i] = sin[\&2] = sin[2\I^3] = sinpvE'J = sm[^2 - 2*3] 

= cos[\E'i + \E'3 ± \E'4] sin[\E'i + \E'3 =F = sin[\&i - \I^2 + ^1^3 ± ^1^4] = (5.28) 







dHsb 


dHsb 






9*3 


9*4 




dHsb 


dHsb 


_ dHsb 


dh 




dh 


dh 



One can therefore fix *j to the discrete values satisfying eqn. ()5.28|) . and only 
solve the remaining 4 equations in (|5.27p . The latter are transformed by the following 
substitutions: 

U4 = y/T4, + K4, M5 = y/T4, - K4, Uq = ^T^ + K5, U-j = y/r^ - K5 (5.29) 

and then multiplied by appropriate factors to remove their denominators. The result 
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is 7 simultaneous polynomial equations: 



dhiui) 

dHsbjuj) 
dhiui) 

dh{ui) 

diAui) 



{uiu^Y'^ = (5.30) 

(MiM2M4M5)^/^ = (5.31) 

{u^u^UeUrY/^ = (5.32) 

{u4U5UQUjy^'^ = (5.33) 



bul + 3ul + 5ul + ^{ul + ul + ul + M7) = P (5.34) 

ul + ul + ul = R (5.35) 

ul-ul + ul~ul = (5.36) 

The 7 unknown variables Ui — u-j are then solved for in PHCpack |l(J2j . 

Treatment of Subsystems When any of the Ui vanishes, at least one of eqns. ()5.3Ut - 
I5.33P becomes unphysical. This happens when there is zero action in any normal 
mode r,-, or when r,- = iK,. 



Special consideration is required for these cases. The contributions from the 
resonance(s) involved with the vanishing action should also vanish from the Hamil- 
tonian. The equations defining the critical points then must be adjusted accordingly. 
As a pedagogical example, in the H2O system the two 0-H normal stretch critical 
point families have no bend action J7j. In deriving these families from the critical 
points analysis, the two stretch-bend Fermi resonances have to be removed. Criti- 
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cal points in the resulting 2 DOF subsystem then correspond to the normal stretch 
modes. 

In order to locate all the possible critical points in the stretch-bend system, all 
combinations of Uj = have to be considered. In each case, the relevant zero-order 
and resonance terms in the classical Hamiltonian are removed. A new canonical 
transformation (Appendix A) may be necessary so that the new angles correspond to 
the remaining resonances. Here we treat the two subsystems outlined in the beginning 
of § 5.3 using this kind of special consideration. 

5.3.2 Results 

The resulting critical points consist of families of curves in the (Tj,Kj) 7- 
dimensional action space, parameterized by the polyad number R. Associated with 
each family is a discrete set of \E'j values. All these solutions turn out to have 
K4 = K5 = 0. In Fig. 15.61 we graph the results in 5 separate panels with each Xj 
versus R separately. In the full system, only those with significantly non-zero Tj are 
displayed here. In the -ft'n/ss subsystem, r2 = r4 = ts = K4 = K5 = 0. In the 
Kn/33 + Ki/2iA subsystem, rs = K4 = = 0. 

With R G [2,8], we found one family for the -ft'n/ss subsystem corresponding 
to the local C-H stretch, two families for the -ft'n/ss + -ft'1/244 subsystem, and two 
families for Hsi, system. 
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FIGURE 5.6 Critical points in the stretch overtone polyads [5R,R,0]. Black, gray 
and red: critical points in if^fe- Navy and magenta: the -ft'n/ss + -^^1/244 system with 
Ts = 0. Green: The i^n/ss system, with T2 = T4 = = 0. 
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When there is only -ft'n/ss resonance, one family of stable (E) critical points 
bifurcates out of the overtone. The new family has 

^1 = (5.37) 

which is the local C-H stretch mode on Fig. 15.31 

When both -R'n/sa + -ft'1/244 are included and ts = K5 = 0, the reduced phase 
space employs four variables (Ji, I2, "^i, '^2)- Two new families of critical points are 
created at: 

(vl/i,v|/2) = (0,0),(0,7r) (5.38) 

The (0, 0) family has (EE) stability and the (0, vr) family has (EH) stability. The C-H 
stretching motion of these critical points is qualitatively like the local C-H stretch 
motion, because the value for both families is the same as in eqn. (j5.37j) . The \l/2 
angle, as the phase angle of a three-mode resonance i^i/244, is yet to be given a clear 
physical meaning. 

In the full Hsh, solving the critical point equations (j5. 30115. 35|) results in two 
more families of critical points with 

(^1, ^2, ^3, ^4) = (0, 0, 0, tt), (0, vr, 0, tt) (5.39) 

Although there is no explicit constraint, both of these families have ^4 = ^5 = 0. The 
(0, 0, 0, vr) family start with (EEEE) stabihty and switches to (EEMM) at = 5.45. 
The (0, TT, 0, tt) family has (EEEH) stability for up to i? = 8. These two families also 



156 

have = 0, indicating a local-type stretch in the C-H DOF. The values of two other 
angles \E'3 = 0, \E'4 = vr, according to Table suggest that the bending motion is of 
the local (L) type. 

A striking feature in Fig. 15.61 is that the critical points in both the -ft'n/ss + 
-^1/244 and full Hamiltonian appear to have bifurcated out of the local C-H stretch 
family (green line). The latter then bifurcates out of the Ts overtones. We believe that 
the new families of critical points in Fig. 15.61 correspond to novel types of dynamics 
in the respective stretch overtone polyads. 

These critical points in the stretch-bend Hamiltonian do not correspond to 
PO in the full phase space when there are up to 3 frequencies in the cyclic angles: 

^'-^ ^^-^^^ 

The last frequency, turns out to be zero for all the critical points graphed in 
Fig. 15. 6( in the same manner as 9^ does in the [Nb-, 0] pure bending polyads when 
£ = (see Chapter 4). The first two frequencies, Op and Or, are not uniquely de- 
fined. Instead, an alternative choice of the canonical transformation can change the 
definitions of and Or as well as their frequencies. 

5.4 Summary 

The critical points analysis is performed on the C2H2 stretch-bend effective 
Hamiltonian, for the purpose of clarifying the fate of normal C-H stretch mode under 
increasing excitation. Preliminary results indicate that the normal antisymmetric C- 
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H stretch overtone (critical point) is first destabilized by Ku/^^ resonance to produce 
the local C-H stretch. Then the local stretch bifurcates into at least 4 families of 
stretch-bend critical points at higher excitation after the inclusion of stretch-bend 
resonances. Such an observation suggests that the stretch-bend dynamics could be 
influences by a chain of bifurcations, each induced by perhaps the addition of one 
resonance. 

Currently we are working to clarify the physical meaning of these critical 
points, as well as their role in the classical phase space. A more complete classification, 
including all the subsystems (combinations of Tj = 0) is also in progress. 
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CHAPTER VI 

CONCLUSIONS AND FUTURE DIRECTIONS 
6.1 Conclusions 

A generalized method of critical points analysis is proposed for studying the 
dynamics of vibrationally excited molecules. The classical form of an effective Hamil- 
tonian with polyad number(s) is canonically transformed to reduce the dimensionality. 
In the reduced phase space, the critical points are systematically found as roots of 
analytic equations. Their number and linear stability are followed as the polyad num- 
ber (s) is varied. Critical points that are linearly stable in all directions are expected 
to indicate regions of regular (quasiperiodic) motion, which correspond to modes of 
vibration. These critical points constitute the most important invariant structure in 
phase space. A change in their number and/or stability, called a bifurcation, indicates 
qualitative changes in both the dynamics and spectra of the system. 

The analysis is carried out in the C2H2 pure bending system. With increasing 
polyad number Nh, four new families of critical points are born in successive bifur- 
cations of the normal mode families. After their first bifurcation, the normal modes 
families become unstable. The two new all-stable families L and CR correspond to the 
new modes dominating the bottom and top (respectively) of the high-lying polyads. 
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The four new families can be qualitatively interpreted as superpositions of resonant 
modes caused when DD-I or £ resonances act separately on the zero-order system. 

The same analysis is extended to the C2H2 stretch-bend system. Preliminary 
results are presented and discussed for polyads [5i?, R, 0] , which contain the stretch 
overtones. As the stretch polyad number R is increased, first the antisymmetric C-H 
normal stretch overtone is first substituted by the local C-H stretch as the stable 
critical point, in a bifurcation induced by i^ii/33 resonance. Then i^i/244 and other 
resonance couplings set in sequentially, creating new families of critical points. These 
critical points are expected to correspond to novel modes of stretch-bend vibration. 

6.2 Summary of Contributions 

In this thesis we have 

• Introduced a generalized critical points analysis method formulated for arbi- 
trary DOF and multiple polyad numbers. It locates new modes of vibration by 
following the critical points in the reduced phase space. The method scales well 
with additional dimensionality of the problem. 

• Found four new modes of vibration in the acetylene pure bending system, and 
provided a qualitative explanation of their origin and nature in terms of single 
DD — I or i resonances. 

• Studied for the first time the acetylene stretch-bend system with all the reso- 
nances. Preliminary results suggest a series of resonances that act sequentially 
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to couple the C-H stretching overtones to the remaining DOF. 

6.3 Future Work 

Chapter 5 only focused on the dynamics of stretch overtone polyads of acety- 
lene. A comprehensive critical points analysis of the full phase space remains to be 
carried out. The resulting critical points could be used to assign the many strongly 
perturbed eigenstates in polyads [4, 20, 0]"^ and [5, 25, 0]"+. These eigenstates ap- 
pear as neither normal mode ZOS, nor attributable to perturbations from a single 
resonance source |117j . 

The method formulated in this thesis has opened the door to understanding the 
dynamics of other high DOF systems with polyad structure. These systems include 
the other C2H2 isotopomers, formaldehyde (H2CO) |125j . methane (CH4) |126j and 
even myoglobin |127j . The last case indicates that even in large biomolecules, under 
favorable conditions a few strongly coupled modes may remain dynamically isolated 
for a (relatively) prolonged time. Analyze these new systems will greatly expand the 
application of the polyad Hamiltonian model. 

New theoretical inquiries also arise from this thesis, especially with regard to 
the mathematical theories of relative equilibria. Most existing studies are concerned 
with rigorously conserved symmetries, such as the angular momentum of an isolated 
body (§ 3.3.2). In the case of molecules, however, both the effective Hamiltonian 
and polyad numbers are approximate. In equating the dynamics of the effective 
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Hamiltonian and that of the true molecular Hamiltonian, the critical points analysis 
needs to be structurally stable with regard to small polyad-breaking terms. These 
terms become increasingly important at high vibrational excitation, especially near 
the threshold of an isomerization barrier |128|129j . A deeper insight into the effect of 
these terms would be of great practical importance to our method. 
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APPENDIX A 
CANONICAL TRANSFORMATION 

This appendix explicitly derives the canonical transformation that results in 
a reduction of the classical Hamiltonian with polyad number (s). The application of 
eqns. fl2.11|) to an A^-mode quantum effective Hamiltonian ()2.1|) results in a classical 
Hamiltonian with N pairs of action-angle variables (tj, (pi) 

H{Ti, (Pi) = Hoin) + H^in, (Pi) (A.l) 

Each resonance coupling in can be expressed as an A^-vector. Let there be 
M (M < A^) resonance vectors that are linearly independent of each other: 

Vi = {Na, N,2, . . . , N,m} for^ = l,...,M (A.2) 

There exist a total number of (A^ — M) polyad numbers j27l28j 

p. = {P,-!, . . . , P,.^} for J = 1, . . . , (AT - M) (A.3) 

which correspond to vectors perpendicular to all Vf. 

Pj ■Vi = (A.4) 

The reduction of eqn. ()A.1|) consists of finding a canonical transformation 
which (1) linearly combines Tj into M new actions Jj and N — M polyad numbers Pj] 



163 

(2) linearly combines 0j into M angles \E'j conjugate to /j and N — M cyclic angle 9j. 
After this transformation, the reduced Hamiltonian will be spanned by only 2N — 2M 
new variables \E'j): 



(A.5) 



This transformation can be expressed in matrix notation as 



9i 



(A.6) 



h [ B 

hi 
Pi 



\ Pn-M ) \ '^N / 

Both A and B are N x N matrices. The first M rows of A are the Vj vectors 



in eqn. ()A.2|) . while the last N — M rows of B are the Pj vectors in eqn. ()A.3|) . 
Using the symplectic formulation of canonical transformations (Chapter 9.3 of 
we require: 



A 
B 



This is equivalent to 



-Em 
En 



A 
B 





En 



-En 




(A.7) 



AB^ = BA^ = E 



N 



(A.8) 
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In solving for the unknown elements in matrices A and B, there are usually 
more unknowns than the number of independent equations. The transformation then 
is not uniquely determined. 

Nevertheless, the classical dynamics (and physical properties in general) should 
not be dependent on the choice of a special coordinate system. Consider a critical 
point in coordinates (/,, where 

dH dH 



(A.9) 



In an alternative coordinate system ( J,, the critical point equations can be derived 



using the chain rule from calculus: 



dH \ / dh 

dJi \ / dJi 

dH j \ dh_ 

dJM / \ dJM 

dH \ / a£i 

a<i>i \ / e<i>i 

dH I \ 



dhi 
dJi 




(A.IO) 



(A.ll) 



d<S>M / \ d<S>M 

The two square matrices are determined by the transformation between (Jj, \l/j) and 
(Jj,$j). From eqn. ()A.9|) . the two column vectors (^fj-j and (^■§§-^ vanish at the 
critical points. The two column vectors on the left, (^fj") (j^) ^^en must also 
vanish. Therefore the critical points are indeed invariant under different choices of 
the canonical transformation. 

Apart from the uncertainty in the canonical transformation, there is also un- 
certainty in choosing the polyad numbers Pj as well as the cyclic angles 6j. The main 
consequence is that the N — M frequencies di associated with a critical point are 



165 

also arbitrary. A preferred definition is not evident from the general consideration ^. 
Additional system-specific constraints, such as relating 9 variables to the vibrational 
frequencies in the Cartesian coordinates, may be necessary to address the choice of a 
non-unique set of action-angle variables. 

In summary, when reducing the classical Hamiltonian using the polyad num- 
ber (s), there is freedom in choosing both the polyad number (s) and the canonical 
transformation defining the reduced phase space. However these different choices 
should lead to the same critical points and physical behavior. 



^The only exceptions are some trivial cases. For example, if a zero-order mode i 
is uncoupled, then using 0j as a cyclic angle is more intuitive. 
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APPENDIX B 

TOPOLOGY OF [Nb, 0] PURE BENDING PHASE SPACE 
B.l. The Poincare-Hopf Index Theorem 

The Poincare-Hopf Index Theorem |13Uj . initially proposed by Poincare and 
later extended by Hopf, provides a constraint on the possible combination of the 
critical points in the reduced classical phase space. 

The theorem states that "The index of a vector field with finitely many zeros on 
a compact, oriented manifold is the same as the Euler characteristic of the manifold" 
PH] . In our cases, the manifold is the reduced phase space while the vector field is the 
flow generated by Hamilton's equations of motion. An index gi is assigned to each 
critical point (i.e. a zero of the vector field) based on its linear stability. The theorem 
indicates that the sum of all indices gi is equal to the Euler characteristics x- X is 
also known as the topological index since it is entirely determined by the topology of 
the manifold. 

A critical point is non-degenerate if its stability matrix A in eqn. ()3.14j) has 
no zero eigenvalues. At such a point, gi = ±1 according to one of the following two 
equivalent criteria: (1) (—1)", with n the number of eigenvalues with positive real 
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part; (2) the sign of the determinant of A. Examples of Qi in 1-3 DOF Hamiltonian 
systems are hsted in Table IB.ll 

TABLE B.l Indices of critical points in Hamiltonian systems with 1-3 DOF. 



Linear Stability 


9i 


Linear Stability 


9i 


E 


+1 


EEE 


+1 


H 


-1 


EEH 


-1 


EE 


+1 


EHH 


+1 


EH 


-1 


HHH 


-1 


HH 


+1 


EMM 


+1 


MM 


+1 


HMM 


-1 



On a 2-dimensional manifold, x is related to the number of "holes" (genus) h 
on the manifold: 

X = 2-2h 

Therefore a torus (one "hole") has x = while a sphere (no "hole") has x = 2. The 
second case has been extensively used by Kellman et al. to verify the critical points on 
a PPS. It proves to be especially useful in determining the stability of critical points 
located at the poles [5D]. Fig. IB . II demonstrates the Poincare-Hopf index theorem, for 
selected PPS in Fig. 13.41 Spheres "1" and "4" each has a total of two stable critical 
points, while sphere "3" has 3 stable and 1 unstable critical points. In all cases, the 
sum of the stability indices by Table equals 2, as expected for a sphere. 
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FIGURE B.l Conservation of topological index on the PPS. A large dot indicates a 
stable (E) critical point and the "X" in sphere 3 indicates an unstable (H) critical 
point. 
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The application of the Poincare-Hopf Theorem to higher dimensional mani- 
folds is limited, mostly due a lack of means of direct visualization. A few known cases 
are summarized in Table lB?2l 

TABLE B.2 Topological indices of 2- and 4-dimensional manifolds, from |131j . 



Symbol 


Manifold 


X 




2- sphere 


2 




2-torus 







4- sphere 


2 




4-torus 


4 



B.2. Topology of [A''^, 0] Bending Phase Space 

This subsection proposes a topological description of the critical points found 
in the C2H2 pure bending polyads [A'j,,0]. 

Currently, there is no sign (such as the display of monodromy |132j ) indicating 
that this classical phase space changes its topology with the Nf,. We first consider 
the low polyad end with < 7, before any bifurcation takes place. In this case 
there are no critical points except possibly where \Ja\ + = Ka. In Figs. l4.lTll4.14l 
of Chapter 4, these locations appear as the boundary of the diamond-shaped space. 
Here this space is schematically illustrated in Fig. IB. 21 Because the coordinate system 
{Ja, Jb, '4'b) IS siugular at these locations, the Hamiltonian is transformed to a local- 
mode representation using the x-K relationship, which is described in Chapter 7.6 of 
[OB]- The result is shown in Fig. IB. 21 It was found that only the four vertices in 
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Fig. IB.2I are critical points: points A and B (J^ = ±Ka, Jb = 0) with (EE) stability 
and points C and D (J^ = 0, = ±Ka) with (MM) stability. 




FIGURE B.2 Critical points in [Nh, 0] polyads at low N^. 
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The two Jfc = families of critical points correspond to states with maxi- 
mum ii in each of the normal mode oscillators ^: 

These two families do not participate in the subsequent bifurcations. 

When Nb > 7, at each bifurcation point in Fig. 14. 5t the sum of indices will 
be conserved if (1) each of the classical normal Trans and Cis motion, being a 2- 
dimensional oscillator, is doubly-degenerate; and (2) each point on the four new 
families of critical points corresponds to a quadruply-degenerate set of PO ^. Then the 
Poincare-Hopf index theorem can be satisfied in the following manner. Substituting 
the stability index of each family and its degeneracy into eqns. fj4.53|) . before and 
after each bifurcation we have: 

L, CR: (EE) ^ 4 (EE) + 2(EH) ; (+2) = 4 x (+1) + 2 x (-1) (B.l) 

Orth: (EH) 4 (EH) + 2(HH) ; (-2) = 4 x (-1) + 2 x (+1) (B.2) 

Pre: (HH) ^ 4 (HH) + 2(EH) ; (+2) = 4 x (+1) + 2 x (-1) (B.3) 

As points C and D do not participate in these bifurcations, only these critical points 
need to be considered. In summary, during each of the four bifurcations of the Trans 
and Cis, the sum of stability indices remains unchanged. 

^Since ii < rii, there cannot be pure "£i -overtones" in the form of e.g. = 0, £4 7^ 

0. 

^In classical mechanics, the two carbon and two hydrogen atoms are assumed 
distinguishable. The quadruply-degenerate PO are related to each other a 90° rotation 
around the C-C bond and a mirror plane perpendicular to it. 
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